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by 


P.  M.  Naghdi 

Department  of  Mechanical  Engineering 
University  of  California,  Berkeley 


Abstract 

The  objective  of  this  paper  is  to  present  an  account  of  recent  develop¬ 
ments  in  the  direct  formulation  of  theories  of  rods  and  shells  based  on 
1  and  2-dimensional  continuum  models  originating  in  the  works  of  Duhem  and 
E.  and  F.  Cosserat.  Following  some  preliminaries  and  description  of 
(3-dimensional)  shell-like  and  rod-like  bodies,  the  rest  of  the  paper  is 
arranged  in  two  parts,  namely  Part  A  (for  shells)  and  Part  B  (for  rods) 
and  can  be  read  independently  of  each  other.  In  each  part,  after  providing 
the  main  ingredients  of  the  direct  model  and  a  statement  of  the  conservation 
laws,  a  rapid  outline  is  given  of  the  derivation  of  the  basic  equations  and 
nonlinear  constitutive  equations  for  elastic  materials.  Each  part  also 
includes  a  discussion  of  constrained  theories  and  an  account  of  recent 
developments  on  the  subject. 
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1.  Introduction 


Rods  and  shells  are  a  class  of  3-dimensional  bodies  whose  boundary 

surfaces  have  special  characteristic  features.  In  general,  two  entirely 

different  approaches  may  be  adopted  for  the  construction  of  1-dimensional 

and  2-dimensional  mechanical  theories  of  rods  and  shells  and  similarly  these 

two  approaches  may  be  used  in  the  construction  of  theories  of  fluid  jets  and 

fluid  sheets.  One  approach  starts  with  the  3-dimensional  equations  of  the 

classical  continuum  mechanics  and  by  applying  approximation  procedures 

strives  to  obtain  1-dimensional  (in  the  case  of  rods  and  jets)  and  2- 

dimensional  (in  the  case  of  shells  and  sheets)  field  equations  and 

constitutive  equations  for  the  medium  under  consideration.  In  the  other 

approach,  the  medium  response  is  modelled  as  a  1-dimensional  and  a  2- 

dimensional  directed  continuum,  called  a  Cosserat  curve  and  a  Cosserat 

surface,  respectively;  and  one  then  proceeds  to  the  development  of  the  field 

* 

equations  and  the  appropriate  constitutive  equations  .  If  full  information 
is  desired  regarding  the  motion  and  deformation  of  the  continuum  under 
study  in  the  context  of  the  classical  3-dimensional  theory,  then  there  would 
be  no  need  to  develop  a  particular  1-dimensional  and  a  2-dimensional  theory. 
In  fact,  the  aim  of  1-dimensional  and  2-dimensional  theories  of  the  type 
mentioned  above  is  to  provide  only  partial  information  in  some  sense:  for 
example,  in  the  case  of  shells  information  concerning  quantities  which  can 
be  regarded  as  representing  the  medium  response  confined  to  a  surface  or 
its  neighborhood  as  a  consequence  of  the  (3-dimensional)  motion  of  the  body, 
or  the  determination  of  certain  weighted  averages  of  quantities  resulting 
from  the  (3-dimensional)  motion  of  the  body. 

_  ™  “ 

Other  2-dimensional  and  1-dimensional  models  may  also  be  used  to  construct 
direct  theories  of  shells  and  rods  but  we  postpone  further  remarks  on 
this  until  later  in  this  section. 
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The  nature  of  the  difficulties  in  the  development  of  both  the  theory 
of  shells  and  the  theory  of  rods  from  the  full  3-dimensional  equations  is 
well  known  and  has  been  elaborated  upon  in  various  contexts  by  Green,  Laws 
and  Naghdi  (1968),  Green  and  Naghdi  (1970),  Naghdi  (1972,  Secs.  1,4;  1974; 
1979a)  and  Ericksen  (1979).  In  view  of  these  it  is  reasonable  to  attempt 
to  formulate  1 -dimensional  and  2-dimensional  theories  of  the  types  described 
above  by  replacing  the  continuum  characterizing  the  (3-dimensional)  medium 
in  question  with  an  alternative  model  which  would  reflect  the  main  features 
of  the  response  of  the  3-dimensional  medium  and  which  would  then  permit  the 
formulation  of  appropriate  1-dimensional  and  2-dimensional  theories  by  a 
direct  approach  and  without  the  appeal  to  special  assumptions  or  approxima¬ 
tions  generally  employed  in  the  derivation  from  the  3-dimensional  equations, 
[t  should  be  emphasized  that  a  Cosserat  surface  and  a  Cosserat  curve  are 
not,  respectively,  just  a  2-dimensional  surface  and  a  1-dimensional  curve; 
but  are,  in  fact,  endowed  with  some  structure  in  the  form  of  additional 
primitive  kinematical  vector  fields. 

The  concept  of  'directed'  or  'oriented'  media  originated  in  the  work 


of  Duhem  (1893)  and  a  first  systematic  development  of  theories  of  oriented 
media  in  one,  two  and  three  dimensions  was  carried  out  by  E.  and  F.  Cosserat 
(1909).  In  their  work,  the  Cosserats  represented  the  orientation  of  each 
point  of  their  continuum  by  a  set  of  mutually  perpendicular  rigid  vectors. 

The  purely  kinematical  aspects  of  oriented  bodies  characterized  by  ordinary 
displacement  and  the  independent  deformation  of  N  deformable  vectors  in 
N-dimensional  space  has  been  discussed  by  Ericksen  and  Truesdell  (1958),  wlu> 
also  introduced  the  terminology  of  directors . 

A  complete  general  theory  of  a  Cosserat  surface  with  a  single  deformable 
director  given  by  Green,  Naghdi  and  Wainwright  (1965)  was  developed  within 
the  framework  of  thermomechanics.  This  derivation  (Green  et  al.  19t>5)  is 
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carried  out  mainly  from  an  appropriate  energy  equation,  together  with 
invariance  requirements  under  superposed  rigid  body  motions.  A  related 
development  utilizing  three  directors  at  each  point  of  the  surface,  in  the 
context  of  a  purely  mechanical  theory  and  with  the  use  of  a  virtual  work 
principle,  is  given  by  Cohen  and  DeSilva  (1966b).  A  further  development 
of  the  basic  theory  of  a  Cosserat  surface  along  with  certain  general 
considerations  regarding  the  construction  of  nonlinear  constitutive  equa¬ 
tions  for  elastic  shells  is  given  by  Naghdi  (1972,  Sec.  8),  which  also 
contains  additional  historical  remarks  relevant  to  oriented  continua  and 
to  the  theory  of  thin  elastic  shells.  Hierarchical  theory  of  Cosserat 
surfaces,  namely  that  comprising  a  material  surface  with  K  (>1)  directors, 
is  contained  in  a  paper  by  Green  and  Naghdi  (1976a)  which  deals  with  fluid 
sheets  and  its  application  to  water  waves. 

A  parallel  development  in  the  theory  of  a  Cosserat  curve  with  two 
deformable  directors  begins  with  a  paper  of  Green  and  Laws  (1966)  whose 
derivation  is  carried  out  mainly  from  an  appropriate  energy  equation, 
together  with  invariance  requirements  under  superposed  rigid  body  motions. 

A  related  development  of  a  directed  curve  with  three  deformable  directors 
at  each  point  of  the  curve,  in  the  context  of  a  purely  mechanical  theory 
and  with  the  use  of  a  virtual  work  principle,  is  given  by  Cohen  (1966).  A 
further  development  of  the  basic  theory  of  a  Cosserat  curve  along  with 
certain  general  developments  regarding  the  construction  of  nonlinear 
constitutive  equations  for  elastic  rods  is  given  by  Green,  Naghdi  and 
Wenner  (1974b).  Hierarchical  theory  of  Cosserat  curves,  namely  that 
comprising  a  material  curve  with  L(>2)  directors,  is  contained  in  a  paper 
by  Naghdi  (1979b)  which  is  concerned  with  applications  to  Newtonian  and 
non-Newtonian  flows  in  pipes. 

Of  course,  the  introduction  of  an  alternative  model  and  formulation  of 
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1-dimensional  and  2-dimensional  theories  by  the  direct  approach  do  not  mean 
that  one  ignores  the  nature  of  the  field  equations  in  the  3-dimensional 
theory.  In  fact,  some  of  the  developments  of  the  field  equations  by  direct 
procedure  are  materially  aided  or  influenced  by  available  information  which 
can  be  obtained  from  the  3-dimensional  theory.  For  example,  the  integrated 
equations  of  motion  from  the  3-dimensional  equations  provide  guidelines 
for  a  statement  of  1-  and  2-dimensional  conservation  laws  in  conjunction 
with  the  1-  and  2-dimensional  models,  and  also  provide  some  insight  into 
the  nature  of  inertia  terms  and  the  kinetic  energy  in  the  direct  formula¬ 
tion  of  the  1-dimensional  and  2-dimensional  theories. 

Inasmuch  as  most  of  the  difficulties  associated  with  the  derivation 
of  the  1-dimensional  and  2-dimensional  theories  from  the  3-dimensional  equa¬ 
tions  occur  in  the  construction  of  the  constitutive  equations,  it  is  in  fact 
here  that  the  direct  approach  offers  a  great  deal  of  appeal.  These  con¬ 
structions,  as  well  as  the  entire  development  by  the  direct  approach,  are 
exact  in  the  sense  that  they  rest  on  (1-dimensional  and  2-dimensional) 
postulates  valid  for  nonlinear  behavior  of  materials  but  clearly  they  cannot 
be  expected  to  represent  all  the  features  that  could  only  be  predicted  by 
the  relevant  full  3-dimensional  equations.  Theories  constructed  via  a 
direct  approach  necessarily  satisfy  the  requirements  of  invariance  under 
superposed  rigid  body  motions  that  arise  from  physical  considerations  and, 
of  course,  they  are  also  consistent  and  fully  invariant  in  the  mathematical 
sense.  Moreover,  the  development  by  the  direct  approach  is  conceptually 
simple  and  does  not  have  the  difficulties  associated  with  approxi mat  ions 
usually  made  in  the  development  of  the  theory  of  thin  shells  or  the  theories 
of  slender  rods  from  their  corresponding  3-dimensional  equations. 

Although  the  direct  approach  to  shells  and  rods  employed  in  this  paper 
is  based  on  the  2-dimensional  and  1-dimensional  directed  continuum  models, 
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respectively,  other  direct  2-dimensional  and  1 -dimensional  models  may  also 
be  used  to  construct  theories  of  shells  and  rods.  For  example  in  the 
case  of  shells,  instead  of  developing  a  theory  based  on  a  Cosserat  surface, 
we  may  consider  only  a  material  surface  and  construct  a  direct  theory  in 
which  the  basic  kinematical  ingredients  are  the  position  vector  of  the 
surface  together  with  its  first  and  second  gradients.  A  theory  of  this 
kind  has  been  discussed  by  Balaban,  Green  and  Naghdi  (1967)  and  a  somewhat 
less  general  theory  by  Cohen  and  DeSilva  (1966a, 1968) .  Although  these 
developments  have  some  overlapping  features  with  corresponding  results  in 
the  theory  of  Cosserat  surfaces,  they  are  more  restrictive.  Additional 
related  remarks  are  made  in  Sec.  6  of  this  paper. 

Following  some  general  background  information  and  definitions  of 
shell-like  and  rod-like  bodies  in  Sec.  2,  the  remained!  of  the  paper  is 
arranged  in  two  parts  which  can  be  read  independently  of  each  other:  one 
part  (Part  A)  is  concerned  with  the  theory  of  shells  and  the  other  (Part  B) 
is  devoted  to  the  theory  of  rods.  In  Part  A  (Secs.  3-8),  first  a  concise 
development  of  the  basic  theory  of  a  Cosserat  surface  with  a  single  director 
followed  by  its  generalization  is  presented.  For  a  Cosserat  surface  with  a 
single  director,  constitutive  equations  are  discussed  in  the  context  of  finite 
deformation  of  elastic  shells  and  a  procedure  is  indicated  for  identification 
of  the  assigned  fields  and  the  inertia  coefficients  which  occur  in  the  basic 
theory.  Next,  a  fairly  detailed  account  of  constrained  theories  of  shells  is 
presented  which  includes  the  construction  of  an  interesting  nonlinear  con¬ 
strained  theory  not  discussed  previously  in  the  literature.  This  is  followed 
by  an  account  of  recent  developments  pertaining  to  elastic  shells  and  a 
representation  of  the  basic  equations  of  a  Cosserat  surface  in  direct 
(coordinate-free)  notation.  A  table  of  contents  for  Part  A  is  listed  in 
the  introductory  paragraph  of  Sec.  3. 
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Similarly,  in  Part  B  (Secs.  9-13),  first  a  concise  development  of  the 
Cosserat  curve  with  two  directors  and  its  generalization  is  presented. 
Next,  with  reference  to  a  Cosserat  curve  with  two  directors,  constitutive 
equations  are  discussed  for  finite  deformation  of  elastic  rods  and  a 
procedure  is  indicated  for  identification  of  the  assigned  fields  and  the 
inertia  coefficients  which  occur  in  the  basic  theory.  This  is  followed 
by  some  additional  remarks  pertaining  to  elastic  rods,  together  with  a 
brief  discussion  of  the  constrained  theories  of  rods,  and  a  representation 
of  the  basic  equations  for  a  Cosserat  curve  in  direct  (coordinate- free) 
notation.  A  table  of  contents  for  Part  B  is  listed  in  the  introductory 
paragraph  of  Sec.  9. 
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2. 


General  Background 

In  this  section,  we  provide  appropriate  definitions  for  shell-like  and 

rod- like  bodies.  To  this  end,  consider  a  finite  three-dimensional  body  £ 

in  a  Euclidean  3-space,  and  let  convected  (or  Lagrangian)  coordinates  01 

(i  =  1,2,3),  be  assigned  to  each  particle  (or  material  po’nt)  of  i.  Further, 

"f*  * 

let  r  be  the  position  vector,  from  a  fixed  origin,  of  a  typical  particle 
of  £  in  the  present  configuration  at  time  t.  Then,  a  motion  of  the  (three- 

dimensional)  body  is  defined  by  a  vector-valued  function  r  which  assigns 

*  '  § 
position  r  to  each  particle  of  £  at  each  instant  of  time,  i.e.  , 

*  ^*123 

r  =  r  (0  ,0  ,0  ,t)  .  (2. 

We  assume  that  the  vector  function  r  --  a  1-parameter  family  of  configura¬ 
tions  with  t  as  the  real  parameter--  is  sufficiently  smooth  in  the  sense  that 
it  is  differentiable  with  respect  to  01  and  t  as  many  times  as  required.  In 

3 

some  developments,  it  is  convenient  to  set  9  =  5  and  adopt  the  notation 

e1  =  (9°,C)  ,  03  =  c  .  (2. 


(2. 

(2. 


We  recall  the  formulas 

/s  * 

3r 


fi 


30 
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T  ’  hj  =  h’h  »  g  =  det(*ij> 


i  ci  i  ii  i  i  ij 

g  *  g*  =  6.  >  g  *  g  g.  .  g  *  g‘  -  g 


dv  =  g^d0  d0  d0‘ 


•4* 

The  use  of  an  asterisk  attached  to  various  symbols  is  for  later  convenience. 
The  corresponding  symbols  without  the  asterisks  are  reserved  for  different 
definitions  or  designations  to  be  introduced  later. 

^ Recall  that  when  the  particles  of  a  continuum  are  referred  to  a  convected 
coordinate  system,  the  numerical  values  of  the  coordinates  associated  with 
each  particle  remain  the  same  for  all  time. 


and  further  assume  that 


g2  =  fglg2i3^  >  0  • 

In  (2.4),  g  and  g1  are  the  covariant  and  the  contravariant  base  vectors  at 
time  t,  respectively,  g„  is  the  metric  tensor,  g1-1  is  its  conjugate,  S* 
is  the  Kronecker  symbol  in  3-space  and  dv  the  volume  element  in  the  present 
configuration. 

★ 

1710  velocity  vector  v  of  a  particle  of  the  three-dimensional  body  in 
the  present  configuration  is  defined  by 


v  =  r 


where  a  superposed  dot  denotes  material  time  differentiation  with  respect 

to  t  holding  0l  fixed.  The  stress  vector  t  across  a  surface  in  the  present 

★ 

configuration  with  outward  unit  normal  v  is  given  by 


1  =  vi  =%  =  V  fk  ’  I  =  gi  *  g~  T  =  T  Si  55  gk  > 


where 


„i  h  ij  h  i  i  *  *  i  *i  i  i  i  j 

T  =gTJg.=gT.g-  ,  v  =  v.g  =  V  g.  ,  rJ=g  •  T  gJ 
J  3"'  v'  t  ~  -'t  — 

ik 

where  T  is  the  symmetric  Cauchy  stress  tensor,  x  its  contravariant  com¬ 
ponents  and®  denotes  the  tensor  product  of  two  vectors.  In  terms  of 
quantities  defined  in  (2. 5) -(2. 8),  the  local  field  equations  which  follow 
from  the  integral  forms  of  the  three-dimensional  conservation  laws  for 
mass,  linear  momentum  and  moment  of  momentum,  respectively,  are 


The  choice  of  positive  sign  in  (2.5)  is  for  definiteness.  Alternatively, 
for  physically  possible  motions  we  only  need  to  assume  that  g?  /  0  with 
the  understanding  that  in  any  given  motion  [ g i g2S3 ]  is  either >0  or <  0 . 
The  condition  (2.5)  also  requires  that  91  be~a~rlght-handed  coordinate 
system. 


(2.5) 


(2.7) 


--V-; 


~5~T2 

P  g  =  0  , 

i  *  *  i,  *  p.  *  :  (2.9) 

t  +  p  f  g'r  =  p  g2v  .  g  *  r  =  o  , 

~  ,  X  ~  '"'A  ~ 

*  * 
where  p  is  the  3-dimensional  mass  density,  f  is  the  body  force  field  per 

unit  mass  and  a  comma  denotes  partial  differentiation  with  respect  to  01 .  For 

later  reference,  we  note  that  for  an  incompressible  medium,  the  condition  of 

incompressibility  may  be  expressed  as 


'2 


g 


★ 

0  or  div  v 


0 


(2.10) 


A  material  surface  in  iS  can  be  defined  by  the  equation  £  =  £(0U);  the 

ct 

equation  resulting  from  (2.1)  with  £  =  £(0  )  represents  the  parametric  form  of 
this  materia]  surface  in  the  current  configuration  and  defines  a  1 -parameter 
family  of  surfaces  in  space,  each  of  which  we  assume  to  be  smooth  and  non¬ 
intersecting.  We  refer  to  the  surface  £=0  in  the  current  configuration  by  s. 

Any  point  of  the  surface  s  is  specified  by  the  position  vector  r,  relative  to 

* 

the  same  fixed  origin  to  which  r  is  referred,  where 

r  =  r(6a,t)  =  r  (ea,0,t)  .  (2.11) 

Ct 

Let  a^  denote  the  base  vectors  along  the  9  -curves  on  the  surface  s.  By 
(2.11)  and  (2.3) 

9r 


a  = 
~a 

and  the  unit  normal  a^  =  a^ 

a  •  a,  =  0  , 

~a  ~3 

We  also  recall  the  formulas 


aotB 

a  aB 

?  =  a 


a  (8^,t)  =  — ~  =  g  (0Y,O,t)  , 

~a  30Ot  ~a  ' 

Y  *  * 

(0  ,t)  to  s  may  be  defined  by 

3 

a3  ’  a3  =  1  ’  ?3  =  a  >  0 

=  a  •  a„  ,  a  =  det(a  „)  , 

-cl  ~8  aB 

a  8  aB  ay  fa 

,  a  •  a  =  a  ,  a  a  „  =  6. 


(2.1 2) 

(2.13) 


(2.14) 


and 


h  .  =  b.  =-a  •  a,  „  =  a  •  a  D  , 
aB  Ba  ~a  ~3,B  -3  .a,B 

Y  ' 

a  1 0  =  b  n  a ,  ,  a,  =  -  b  a  ,  b  _ i  =  b  i _  , 

-a[B  aB-3  ~3,«  a~y  aB|y  ay|B 

The  use  of  the  same  symbols  for  base  vectors  of  a  surface  in  (2 .  1 2) -  (2 . 1 3) 
and  for  the  triad  of  a  space  curve  in  (2 . 1 7) - (2 . 1 8 )  should  not  give  rise  to 
confusion.  The  main  developments  for  shells  and  rods  arc  dealt  with 
separately  in  the  rest  of  the  paper;  this  permits  the  use  of  the  same  symbol 
for  different  quantities  in  the  case  of  shells  and  rods  without  confusion. 


where  a  denote  the  reciprocal  base  vectors  of  the  surface  s,  a  D  and  b  arc 

~  Otp  Ctp 

the  components  of  its  first  and  second  fundamental  forms,  a  comma  denotes  partial 
differentiation  with  respect  to  the  surface  coordinates  0^,  a  vertical  bar  stands 


for  covariant  differentiation  with  respect  to  a^  and  6^  is  the  Kronecker  symbol 
in  2- space. 

A  material  line  (not  necessarily  a  straight  line)  in  $  can  be  defined  by 
the  equations  0a  =  0a(£);  the  equation  resulting  from  (2.1)  with  0a  =  0a(£)  repre¬ 
sents  the  parametric  form  of  this  material  line  in  the  current  configuration 
and  defines  a  1-parameter  family  of  curves  in  space,  each  of  which  we  assume  to 
he  smooth  and  nonintersecting.  We  refer  to  the  space  curve  0  =  0  in  the  current 
configuration  by  c.  Any  point  of  this  curve  is  specified  by  the  position  vector 


r,  relative  to  the  same  fixed  origin  to  which  r  is  referred,  where 

r  =  r(£,t)  =  r  (0,0,£,t)  . 

§ 

Let  a^  denote  the  tangent  vector  along^the  £-curve.  By  (2.16)  and  (2.3)^, 

3r 

a3  =  a3(£,t)  =  ^  =  g3(0,0,C,t) 

and  the  uni,.  principal  normal  a^  and  the  unit  binormal  vector  a 2  to  c  may  be 


(2.16) 


(2.17) 


introduced  as 


?1  = 


9a37H 


,  a2  =  =  ~|~a  |  *  al  ’ 


(2. IK) 


I  hi  “  (a33)Z  ’  a33  =  ?3>h  *  [?ia2?3]  >  °  * 

where  the  notation  ja3|  stands  for  the  magnitude  of  a^.  The  system  of  base 

vectors  a^  are  oriented  along  the  Serret-Frenet  triad  and  satisfy  the  dif¬ 
ferential  equations 

3a  i  j,  9^2  1  ^a33 

rtf;'  =  T(a33)  'a2  *  Ka3  ’  TT  =  ’  Tl'a33)  ’  “9F  =  a33K?l  +  2a^  K  ~3 


(2.19) 


where  <  and  t  denote,  respectively,  the  curvature  and  the  torsion  of  c.  !n  the 
special  case  that  c  is  a  plane  curve,  we  may  choose  a^  as  the  unit  normal  to  the 
curve  and  then  a^  will  be  perpendicular  to  the  plane  of  and  ay  If  c  is  a 
straight  curve,  then  there  is  no  unique  Serret-Frenet  triad  and  a^  may  be  chosen 
as  any  orthogonal  triad  with  as  unit  vectors.  Equations  (2.19)  are  not 

^The  designation  of  the  tangent  vector  to  a  curve  by  a3  should  not  be  confused 
with  the  use  of  the  same  symbol  for  a  different  purpose  in  (2.13).  In  this 
connection,  see  the  preceding  footnote. 
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identical  to  the  formulas  of  Frenet  because  the  parameter  F  is  not  necessarily 
the  arc  length  of  c.  It  may  be  noted  here  that  the  convected  coordinate  f,  may 
be  chosen  to  coincide  with  the  arc  length  in  any  one  configuration  of  the 
material  curve,  e.g.,  in  the  present  configuration.  However,  in  a  general 
motion  (involving  different  configurations)  the  arc  length  between  any  pair  of 
particles  changes  while  the  convected  coordinates  of  each  particle  must  remain 
the  same.  Therefore,  arc  length  would  not  qualify  as  a  convected  coordinate. 

In  the  next  four  paragraphs  (identified  as  subsections  2A  and  2B)  we 
provide  appropriate  definitions  for  shell-like  and  rod-like  bodies  in  fairly 
precise  terms. 

2A.  Definition  of  a  shell-like  body.  A  representation  for  the  motion  of 

a  thin  shell. 

Consider  a  two-dimensional  surface  s  defined  by  the  parametric  equa - 
tion  F.  =  0,  over  a  finite  coordinate  patch  a/<0  <a”,  <g".  Let  r 

and  a^  denote,  respectively,  the  position  vector  and  the  unit  normal  to  s. 

At  each  point  of  s,  imagine  material  filaments  projecting  normally  above 
arid  below  the  surface  s.  The  surface  formed  by  the  material  filaments  con¬ 
structed  at  the  points  of  the  closed  boundary  curve  of  s  is  called  the  lateral 
surface.  Such  a  3-dimensional  body  (depicted  in  Fig.  1)  is  called  a  shell  if  the 
dimension  of  the  body  along  the  normals,  called  the  height  and  denoted  by 
h,  is  smal 1 .  A  shell  is  said  to  be  thin  if  its  thickness  is  much  smaller 
than  a  certain  characteristic  length  L(s)  of  the  surface  s,  for  example,  the 
local  minimum  radius  of  curvature  of  the  surface,  or  the  smallest  dimension 
of  s  in  the  case  of  a  plane.  If  h  is  constant,  the  shell  is  said  to  be  of 
uniform  thickness,  otherwise  of  variable  thickness.  Since  a  material  surface 
in  tiie  three-dimensional  body  can  be  defined  by  the  equation  F  =  FiO  ),  it 
follows  that  the  equation  resulting  from  (2.1)  with  F  =  F(0a)  represents  the 
parametric  form  of  the  material  surface  in  the  present  configuration.  In 
particular,  the  equation  F  =  0  defines  a  surface  in  space  at  time  t,  which 
we  assume  to  be  smooth  and  nonintersecting.  livery  point  of  this  surface 

1  1 . 


I 


I 

I 

4 


has  a  position  vector  r  specified  by  (2.11).  Let  the  boundary  of 
the  three-dimensional  continuum  be  specified  by  the  material 
surfaces 

£  =  ce1,e2)  ,  5  =  c2(e1,e2)  ,  <  c2  ,  (2.2oj 

with  the  surface  £ =  0  lying  either  on  one  of  the  two  surfaces  (2.20),  7  or 
between  them  (see,  for  example,  Fig.  1),  and  a  material  surface 

ffe^e2}  =  o  ,  (2.211 

which  is  chosen  such  that  £  =  const,  form  closed  smooth  curves  on  the  surface 
(2.21).  As  pointed  out  previously  by  Naghdi  (1975a),  in  the  development  of  a 
general  theory,  it  is  preferable  to  leave  unspecified  the  choice  of  the  rela¬ 
tion  of  the  surface  s  (£ =  0)  to  the  major  surfaces  s+  and  s'.  In  special 
cases  of  the  general  theory  or  in  specific  applications,  however,  it  is 

necessary  to  fix  the  relation  of  s  to  the  surfaces  (2.20)  _. 

1  ,  / 

A  * 

We  now  suppose  that  r  in  (2.1)  can  be  represented  by  the  Taylor 
expansion  in  the  bounded  region  £^<£<£2  with  coefficients  which  are 
continuous  functions  of  8a,t  and  have  continuous  space  and  time  derivatives 
of  order  2.  Thus,  for  shell-like  bodies,  we  write 


r*  =  r  +  Z  CNd  ,  d  =  d  (6a,t) 
N=1 


(2.22] 


and  by  (2.3)^  and  (2.6)  we  also  have 


K  N  9d  K 

_  _N  ~N  _  ,,,-N-l. 

~a  =  ~a +  ^  ^  a  ’  §3  =  ^  ^  ’ 

a  “  N=1  80a  0  N=1  ~N 


(2.23) 
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•  T-  - 


JldL 


r 


K 


-N 


Z  =  Z*  Z  ^  WN  '  *  =  r  ’  WN  =  dN  ’ 

N=1  ~ 

where  r  is  defined  by  (2.1)  and  a  superposed  dot  in  (2.22)  denotes  material 
time  differentiation  with  respect  to  t  holding  8  fixed.  A  special  case 
of  (2.22)  which  is  of  particular  interest  in  subsequent  developments  is  when 
N  =  1 ,  namely 


r  +  £d  , 


( 


where  we  have  set  dj  =  d. 


2B.  Definition  of  a  rod-like  body.  A  representation  for  the  motion  of 
a  slender  rod. 

OL 

Consider  a  space  curve  c  defined  by  the  parametric  equations  0  =0, 
over  a  finite  interval  ^<£<£2*  Let  r  be  the  position  vector  of  any  point 
of  c  and  let  a^,a2  and  a^  denote  its  unit  principal  normal,  unit  binormal  and 
the  tangent  vector,  respectively.  At  each  point  of  c,  imagine  material  fila¬ 
ments  lying  in  the  normal  plane,  i.e.,  the  plane  perpendicular  to  a^,  and 
forming  the  normal  cross-section  The  surface  swept  out  by  the  closed 

boundary  :urve  dJl  of is  called  the  lateral  surface.  Such  a  3-dimensional 
n  n 

body  (depicted  in  Fig.  2)  is  called  rod-like  if  the  dimensions  in  the  plane  of 
normal  cross-section  are  small  compared  to  some  characteristic  dimension 
L(c)  of  c  (see  Fig.  2),  e.g.,  its  local  radius  of  curvature  1/k,  or  the 
length  of  c  in  the  case  of  a  straight  curve.  A  rod-like  body  is  said  to 
be  slender  if  the  largest  dimension  of^n  is  much  smaller  than  1(c).  If 
.4  is  independent  of  £,  the  body  is  said  to  be  of  uniform  cross-section, 
otherwise  of  variable  cross-section.  Since  a  material  curve  in  the  three- 
dimensional  body  $  can  be  defined  by  the  equations  0a=Oa(£),  it  follows 
that  the  equation  resulting  from  (2.1)  with  0a  =  0a(C)  represents  the 
parametric  form  of  the  material  curve  in  the  present  configuration  and 
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defines  a  curve  c  in  space  at  time  t,  which  we  assume  to  be  sufficiently 
smooth  and  nonintersecting.  Every  point  of  this  curve  has  a  position 
vector  specified  by  (2.14).  Let  the  (3-dimensional)  rod-like  body  in 
some  neighborhood  of  c  be  bounded  by  material  surfaces  C  = 

(indicated  in  Fig.  2)  and  a  material  surface  of  the  form 

Ffe1,©2^)  =  0  ,  (2.26) 


which  is  chosen  such  that  £  =  constant  are  curved  sections  of  the  body 
bounded  by  closed  curves  on  this  surface  with  c  lying  on  or  within  (2.26) 
In  the  development  of  a  general  theory,  it  is  preferable  to  leave  unspeci¬ 
fied  the  choice  of  the  relation  of  the  curve  c  to  one  on  the  boundary 
surface  (2.26).  In  special  cases  or  in  specific  applications,  however,  it 

is  necessary  to  fix  the  relation  of  c  to  the  surface  (2.26). 

~  * 

We  now  suppose  that  r  in  (2.1)  can  be  represented  by  the  Taylor 
expansion  in  the  bounded  region  lying  inside  the  surface  (2.26)  and 
between  £  =  C  =  C2>  with  coefficients  which  are  continuous  functions 
of  £,t  and  have  continuous  space  and  time  derivatives  of  order  2.  Thus, 
for  rod- like  bodies,  we  write 


r+  E  6 
~  N=1 


aN 

•9  ~a. 


(S,t) 


and  by  (2 . 3) ^  and  (2.6)  we  also  have 


K  ou  cr, 

g  =  d  +  E  N0  . . .9  % 

-8  ~8  N=1  ~Ba2...aN 


g,  =  a,  +  E  9  . 

N=1 


.9  (H 


/3f.)  ,  (2.28) 
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*  K  cx. 

v=v+£8...0w  ,  w  =d  , 

~  ~  N=1  'V*-aN  'V‘-aN  ^r--aN 

where  r  in  (2.27)  is  defined  by  (2.16),  d  is  symmetric  with  respect 

1"  N 

to  indices  a^...a  and  a  superposed  dot  in  (2.29)  denotes  material  time 
differentiation  with  respect  to  t  holding  £  fixed.  A  special  case  of 
(2.27)  which  is  of  particular  interest  in  subsequent  developments  is  when 
N  =  1 ,  name 1 y 


(2 . 2‘.) ) 


r  =  r  ♦  0  d 


(2.. Vi  | 


where  we  have  put  =a. 


*w/>.  X.-W'mt*  ' 


Part  A 


Elastic  shells:  A  direct  formulation 

In  Part  A  (Secs.  3-8),  we  summarize  the  main  kinematics  and  the  basic 
principles  of  the  theory  of  Cosserat  (or  directed)  surfaces  and  then  discuss 
the  constitutive  equations  for  elastic  shells,  as  well  as  several  related 
aspects  of  the  basic  theory  and  recent  developments  on  the  subject. 

Although  we  are  concerned  here  mainly  with  the  purely  mechanical  theory 
involving  appropriate  forms  of  the  conservation  laws  for  mass,  linear 
momentum,  director  momentum  and  moment  of  momentum,  we  also  include  a  state¬ 
ment  of  the  conservation  of  energy.  The  latter  provides  motivation  in  the 
development  of  certain  constitutive  equations,  such  as  those  for  an 
elastic  material,  and  in  the  discussion  of  aspects  of  some  special  solu¬ 
tions  involving  jump  in  energy.  The  contents  of  Part  A  are  as  follows: 

3.  The  basic  theory  of  Cosserat  surfaces 

3.1  Kinematics  of  a  Cosserat  surface  C . 

3.2  Basic  principles  of  a  Cosserat  surface  C. 

3.3  Hierarchical  theories  of  Cosserat  surfaces. 

4.  Elastic  shells. 

5.  Identification  of  the  assigned  fields  and  the  inertia  coefficients. 

6.  Constrained  theories  of  shells 

6.1  Incompressible  Cosserat  surface  C. 

6.2  A  constrained  theory  with  director  along  the  normal 
to  the  surface  of  C. 

7.  Additional  remarks  on  shells. 

8.  Basic  equations  for  a  Cosserat  surface  in  direct  notation. 
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3.  The  basic  theory  of  Cosserat  surfaces 

Having  introduced  the  notion  of  a  (three-dimensional)  shell-like  body 
in  section  2,  we  now  formally  define  a  direct  model  for  such  a  body.  Thus, 
deformable  media  which  are  modelled  by  a  material  surface  S  embedded  in  a 
Euclidean  3-space,  together  with  K  (K = 1,2, . . . ,N)  deformable  vector  fields  -- 
called  directors  --  attached  to  every  point  of  the  material  surface  are  called 
Cosserat  surfaces  or  directed  sur faces  and  may  be  conveniently  referred  to 
as  f  .  The  directors  which  are  not  necessarily  along  the  unit  normals  to  the 

iv 

surface  have,  in  particular,  the  property  that  they  remain  unaltered  in 

length  under  superposed  rigid  body  motions. 

In  the  absence  of  the  directors,  we  merely  have  a  2-dimensional 

material  surface  S  which  can  serve  as  a  model  for  the  construction  by  direct 

approach  of  the  membrane  theory  of  shells.  With  K =  1 ,  the  directed  medium 

is  a  body  -  C  comprising  a  material  surface  and  a  single  deformable  director 

attached  to  every  point  of  the  material  surface  of  C.  The  latter  is  the 

simplest  model  for  the  construction  of  a  general  bending  theory  of  thin 

shells;  and,  for  simplicity,  we  restrict  attention  to  this  particular  model 

* 

in  most  of  the  development  of  section  3. 

3. 1  Kinematics  of  a  Cosserat  surface  C. 

Let  the  particles  of  the  material  surface  S  of  C  be  identified  by  means 

of  a  system  of  convected  coordinates  0a  (a  =1,2)  and  let  the  2-dimensional 

region  occupied  by  the  material  surfaced  in  the  present  configuration  of 

at  time  t  be  referred  to  as  4.  Let  r  and  d  denote  the  position  vector  of  a 

typical  point  of  4  and  the  director  at  the  same  point,  respectively.  Also, 

ot 

let  aa*a^  designate,  respectively,  the  base  vectors  along  the  0  -curves  on  4 

*  * 

A  brief  account  of  the  more  general  theory  for  Cosserat  surfaces  is 

indicated  at  the  end  of  this  section. 
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and  the  outward  unit  normal  to  4.  Then,  a  motion  of  the  Cosserat  surface  is 
defined  by  vector-valued  functions  which  assign  position  r  and  director  d  to 

★  Hf 

each  particle  of  C  at  each  instant  of  time,  i.e.  , 


where 


r  =  r(0  ,t) 


d  =  d(0  ,  t) 


>  0 


a 

~a 


a  (0a,t)  = 
~a 


8r 


(3.1) 


(3.2) 


and  the  condition  (3.1)^  ensures  that  the  director  d  is  nowhere  tangent  to 
4.  The  velocity  and  the  director  velocity  vectors  are  defined  by 


v 


w  =  d 


(3.3) 


and  since  the  coordinate  curves  on  4  are  convected  from  (3.2),  we  have 

a  =  v  , 

-a  ~  ,a 

where  a  superposed  dot  denotes  differentiation  with  respect  to  t  holding 
9™  fixed. 

It  is  convenient  to  introduce  here  a  slightly  different  notation  than 
that  adopted  in  Naghdi  (1972)  and  a  number  of  earlier  papers  on  the  subject. 
Thus,  we  put 


d  =  a  ,  d.  =  d 
-a  ~a  3 

and  observe  that,  in  view  of  (3.1)^  and  (3.4),  are  linearly 

independent  vectors.  Hence,  we  may  introduce  a  set  of  reciprocal  vectors 


(3..S) 


*  * 

For  convenience,  we  adopt  the  notation  for  r  in  (2.11)  and  (2.25)  also  for 
the  surface  (3.1)j.  This  permits  an  easy  identification  of  the  two  surfaces, 
if  desired.  The  choice  of  positive  sign  in  (3.1)j  is  for  definiteness. 
Alternatively,  it  will  suffice  to  assume  that  [aj^dJ^O  with  the  under¬ 
standing  that  in  any  given  motion  the  scalar  triple  product  [; | a  2d ]  is 
either  >0  or  <  0. 
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(3.0) 


d1 


such  that 


d. 

~i 


dj 


where  6?  is  the  Kronecker  symbol  in  3-space.  Whenever  desirable,  the 

notations  d^  =  (d^.d^.d^)  and  (a^.d)  will  be  used  interchangeably  throughout 

Part  A  depending  on  the  particular  context.  Consider  now  a  reference 

configuration,  not  necessarily  the  initial  configuration,  of  the  Cosserat 

surface  C.  In  the  reference  configuration,  let  the  material  surface  of  C 

be  referred  to  by  £  with  R  as  its  position  vector;  let  D  be  the  director  at 

Ct 

R;  and  let  A^,A^  denote,  respectively,  the  base  vectors  along  the  0  -curves 

on  .5  and  the  unit  normal  to  S  .  Then,  in  the  reference  configuration  we  have 
K  K 


R  =  R(0a)  ,  D  =  D(0a)  ,  [A  A2D]  >  0 


(3.7) 


where 


6a  *  6a‘0Y>  -  dr 

99 

and  (3.7)3  ensures  that  D  is  nowhere  tangent  to  the  surface  d>R.  If  the 
reference  configuration  of  C  is  specified  to  be  the  initial  configuration, 
say  at  time  t=0,  then  the  vector-valued  functions  on  the  right-hand  sides 

A  /\ 

Uf  (3.7)  9  can  be  identified  with  r(0a,O)  and  d(0a,O),  respectively. 

*  »  *-  ~  " 

Analogously  to  (3.5),  we  set 


9R 


and  note  that 


D  =  A  ,  D. 
~a  ~a  ~3 


the  dual  of  (3.6)  is  given  by 


D 


(3.8) 


(3.9) 


(3.10) 
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5.2  Basic  principles  of  a  Cosserat  surface  C. 


In  the  development  of  this  subsection,  we  follow  the  mode  of  derivation 
of  the  basic  theory  of  a  Cosserat  surface  employed  by  Naghdi  (1972, 

Sec.  8).  Let  <?,  bounded  by  a  closed  curve  3£P  ,  be  a  part  of  4  occupied 
by  an  arbitrary  material  region  of  S  in  the  present  configuration  at  time 
t  and  let 


v 


a 

v  a~ 
~a 


=  v  a 
cu 


(5.11  j 


be  the  outward  unit  normal  to  3£?  .  It  is  convenient  at  this  point  to  define 
certain  additional  quantities  as  follows:  The  mass  density  p  =  p  (0  ,t)  of 
the  surface  4  in  the  present  configuration;  the  contact  force  n  =  n(9',t  ;  v) 
and  the  contact  director  force  m=m(0  ,  t  ;  v) ,  each  per  unit  length  of  a 
curve  in  the  present  configuration;  the  assigned  force  f=f(0Y,t)  and  the 
assigned  director  force  S,  =  JL(0^,t) ,  each  per  unit  mass  of  the  surface  4; 
the  intrinsic  director  force  k  per  unit  area  of  4;  the  inertia  coefficients 

1  i  y  2  2  Y 

v  =y  (O')  and  y  =  y  (01)  which  are  independent  of  time;  the  specific  internal 
energy  £=e(0Y,t);  the  heat  flux  h  =  h(0Y,t;v)  per  unit  time  and  per  unit 
length  of  a  curve  39*;  the  specific  heat  supply  r=r(0^,t)  per  unit  time;  and 
the  element  of  area  do  of  the  surface  9\  and  the  line  element  ds  of  the  curve  39*. 


The  assigned  field  f  may  be  regarded  as  representing  the  combined  effect  of  (i) 
the  stress  vector  on  the  major  surfaces  of  the  shell-like  body  denoted  by  f  , 
e.g.,  that  due  to  the  ambient  pressure  of  the  surrounding  medium,  and  (ii)  an 
integrated  contribution  arising  from  the  three-dimensional  body  force  denoted 
by  f^,  e.g.,  that  due  to  gravity.  A  parallel  statement  holds  for  the  assigned 


The  notations  for  the  contact  force  n,  the  contact  director  force  m  and 
the  surface  director  force  k  are  the  same  as  those  in  Naghdi  (1977),  but 
differ  from  Naghdi  (1972)  and  most  of  the  previous  papeTs  on  the  subject. 

In  fact,  the  vector  fields  n,m,k  of  Part  A  of  the  present  paper  correspond, 
respectively,  to  N,M,m  in  Naghdi  (1972)  and  most  of  the  previous  papers  on 
the  subject.  Also  the  notations  for  the  inertia  coefficients  y*  and  y^, 
which  occur  in  (3. 1 3) - (3 . 14) ,  differ  from  the  corresponding  notations  in 
previous  papers.  In  most  of  the  previous  papers  (for  example,  Green  and 
Naghdi  1976a,  Naghdi  1975a  or  Naghdi  1979a)  the  notations  M,k“  or  aij  ,«2 
were  used  in  place  of  y^,y2. 
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.A 


field  l.  Similarly,  the  assigned  heat  supply  r  may  be  regarded  as  representing 
the  combined  effect  of  (i)  heat  supply  entering  the  major  surfaces  of  the 
shell- like  body  from  the  surrounding  environment,  denoted  by  r  ,  and  (ii)  a 
contribution  arising  from  the  three-dimensional  heat  supply,  denoted  by  r^. 
Thus,  we  may  write 


f  =  f,  +  f  ,  fc  =  H.  +  Z  ,  r  =  r  +  r 
~b  ~c  ~  ~b  ~c  be 


We  assume  that  the  kinetic  energy  of  the  Cosserat  surface  C  per  unit 


(3.12) 


area  of  4  in  the  present  configuration  is  given  by 


1  2 

«  =  %)(v,v  +  2yv*w  +  yw,w) 


We  further  define  the  momentum  corresponding  to  the  velocity  v  and  the 
director  momentum  corresponding  to  the  director  velocity  w  by 

3k  ,  1,  3k  ,12. 

97  =  p(y*y«)  .  ^  =  p(y  v  +  y  w)  • 


Also,  the  physical  dimensions  of  p,n,f  are 


phys.  dim.  p  =  [ML  ]  , 


phys.  dim.  n  =  [MT  ,  phys.  dim.  f  =  [LT  , 


(3.  1.3) 


(3.14 


(3. 15) 


where  the  symbols  [L],  [M]  and  [T]  stand  for  the  physical  dimensions  of  length, 

mass  and  time.  The  dimensions  of  the  vector  fields  m,£  and  k  depend  upon  the 

*  *  . 

physical  dimension  of  d.  Here  we  choose  d  to  have  the  dimension  of  length 

and  then  m,5L  will  have  the  same  physical  dimensions  as  n,f  in  (.3.15)  while  k 

-1-2 

will  have  the  physical  dimension  of  [ML  T  ]. 


Depending  on  the  choice  of  the  physical  dimension  of  d  and  with  reference 
to  m , l  and  k  the  terminologies  of  the  contact  director  couple,  the  assigned 
director  couple  and  the  intrinsic  director  couple,  respectively,  are  also 
used  in  the  literature.  In  particular,  the  latter  terminologies  are 
employed  in  Naghdi  (1972),  where  d  is  taken  to  be  dimensionless. 
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In  terms  of  the  above  definitions  of  the  various  field  quantities  and 

with  reference  to  the  present  configuration,  the  conservation  laws  in  the 

* 

purely  mechanical  theory  of  a  Cosserat  surface  C  are 


jr  /  p  [c+k  ]do  =  /  p(f  "  v  +  &  *  w  +  r)do  +  /  (n  •  v  +  m  •  w  -  h)ds  .  (  3. 1  7) 

9  9  ~  ~  ~  ~  ~  ~  ~  ~  ~ 

The  basic  structure  of  (3.16)  _  and  (5.17)  and  their  forms  are 

analogous  to  the  corresponding  conservation  laws  of  the  classical  3-dimensional 

continuum  field  theory.  The  structures  of  (3.16)^  and  (3.16)^  are  less 

obvious,  but  a  motivation  for  their  forms  is  provided  by  a  derivation  of 

the  basic  field  equations  for  shell-like  bodies  obtained  from  the 

3-dimensional  equations  of  continuum  mechanics  in  which  the  position  vector 
★ 

r  in  3-space  is  approximated  by  an  expression  of  the  form  (2.19).  It 
should  be  noted  here  that  the  conservation  laws  (3. 1 6) - (3. 1 7)  are  consistent 
with  the  invariance  conditions  under  superposed  rigid  body  motions,  which 


As  the  integrals  on  the  left-hand  sides  of  (3.16)-,  ^  ^  allow  for  coupling 
in  inertia  terms,  they  are  slightly  more  general  than  the  corresponding 
expressions  in  Naghdi  (1972).  The  conservation  laws  (3.16)  with  coefficients 
y*  =  0  and  y-  -  'i  l  0  reduce  to  those  given  by  F.qs.  (8.111  in  Naghdi  (1972). 


T  ■> 


ordinarily  have  wide  acceptance  in  continuum  mechanics.  Moreover,  as  shown 


in  Naghdi  (1972,  Sec.  8),  the  conservation  laws  (3.16)  ,  (3.16)2  and  ( 3 . 1  C> )  are 
equivalent  to,  and  can  be  derived  from  the  conservation  of  energy  (3.17) 
and  the  invariance  conditions  under  superposed  rigid  body  motions.  The 
conservation  law  (3. lb)  for  the  director  momentum  must  be  postulated 
separately. 

Returning  to  the  conservation  laws  (3.16)  and  (3.17),  we  note  that 
under  suitable  continuity  assumptions  the  contact  force  n,  the  contact 
director  force  m  and  the  heat  flux  h  can  be  expressed  in  the  forms  (for 
details  see  Naghdi  1972,  Sec.  8): 


..a 

n  =  N  v 


a 


m 


=  iA 


a 


a 

q  v 


Ct  n  01 

where  N  ,t/r  transform  as  contravariant  surface  vectors  and  q  are  the 
contravariant  components  of  the  heat  flux  vector 


q  =  q  a 


a 


With  the  use  of  (3.18)  and  by  usual  procedures,  from  the  conservation  laws 
(3.15)  and  (3.16)  follow  the  local  field  equations 


and 


who  re 


*  Ct 

pa 2  =  A  or  p  +  pa  •  v  Q  =  0  > 


Na,  +  pf  =  p(v  +y1w)  , 


M01 1  +  pS,  -  k  =  pty^v+y^w)  , 


a  x  Na  +  d  k  k  *  d  x  M01  =0 
a  ~  -  ~,a  ~ 


pr  -  qa,  -  pc  +  P  =  0  , 


(3.18) 


(3.  I  9 ) 


(3.20) 


(3.31  ) 


23. 


(3.22) 


P 


v  + 

~  .a 


k  •  w  +  M  •  w 


Ct 

is  the  mechanical  power,  X  in  (3.19)^  is  a  function  of  0  only,  a  comma 

Ct 

denotes  partial  differentiation  with  respect  to  0  ,  a  vertical  line  stands 
for  covariant  differentiation  with  respect  to  the  metric  tensor  of  the 
surface  4  and 


^1-2-3^ 


(3.23) 


3. 3  Hierarchical  theories  of  Cosserat  surfaces 

Although  the  theory  outlined  in  subsection  3.2  is  sufficiently  general 
for  many  applications,  on  occasion  it  becomes  necessary  to  consider 
Cosserat  surfaces  with  more  than  one  director.  Therefore,  we  now  briefly 
discuss  the  kinematics  and  the  balance  laws  of  Cosserat  surfaces  C 

K 

having  K  (K=l,2,...)  directors  attached  to  every  point  of  a  material 
surface  S.  Thus,  we  admit  K  directors  at  r  denoted  by  dM  (M =  1 , 2, . . . , K) ; 
and,  instead  of  (3.1).  _  specify  a  motion  of  ('  by 

1  j 

r  =  r(0a,t)  ,  dM  =  dM(6a,t)  >. 

The  velocity  vector  is  still  given  by  (3. 3) ^  but  corresponding  to  (3.3)7 
we  now  define  the  director  velocities 

~M  = 

We  recall  for  K=1  (C^-O,  the  kinetical  quantities  introduced 


(3.24) 
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in  subsection  3.2  consisted  of  n,k,m  and  the  assigned  fields  f ,£.  Keeping 
this  in  mind,  for  a  body  C ^  we  admit  more  general  kinetical  quantities  and 
assigned  fields 


.  N  N 
n  ,  k  ,  m  , 


f  ,  A 


( 3 . 2(i 


for  N  =  1 ,2, . . . ,K,  and  corresponding  to  (3.13)  and  (3.14)  write  the  more 

l  9 

general  expressions  for  kinetic  energy  of  C  and  associated  momentum  and 
director  momentum,  namely 


<  =  ho 


£  M+N 

l  y  w. 

M,  N=o 


w 

M  ~N 


=  v 


3k 

3v 


p(v  + 


;  m  . 

1  y  = 

M=1 


V  M 
p  l  y  w 

M=o 


M 


(3.271 


_3k 

3w, 


=  p(y  v  + 


N 


K 

Z 

W*1 


M+N  . 

y  ~m^  = 


K 

p  Z  y 
M=o 


M+N 


each  per  unit  area  of  the  surface  2*.  The  inertia  coefficients  yM+N  are 

Ct 

functions  of  0  only  and  satisfy  the  conditions 


M+N  N+M  M+o  o+M  M  o 

y  =  y  ,  y  =  y  =  y  ,  y  =  i  •  (3.28) 


In  the  special  case  of  <?  (=  C)  we  may  use  the  notations 

2l  =  $  '  "l  =  «  •  (S--0) 

For  a  detailed  statement  of  conservation  laws  appropriate  for  Cosserat 

surfaces  we  refer  the  reader  to  Green  and  Naghdi  (1976a,  Sec.  2)  but 

indicate  here  the  structure  of  the  corresponding  local  field  equations.  In 

this  connection,  we  first  note  that  for  a  purely  mechanical  theory  by  usual 

procedure  in  addition  to  (3.18).  we  now  obtain  mN =  MN°‘v  .  Then,  the  local 

J  a 

field  equations  for  Cosserat  surface  are: 
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Ay  =  o  ,  A  =  o  ,  A  =  pa*  , 


..a  .  1  M  • 

"  |a  +  p!  *  p.tZ  y  ~M  ’ 
'  M=o 


Jia  oN  v  M+N  '  .  - 

M  ia  +  p£  -  k  =  p  Z  y  w  ,  (N  =  1 , 2 , .  . . ,  K)  , 

'  ~  M=o 

K  K 

a/N“t  E  d..xk%  Z  d..  x  M™*  =  o  . 

~a  ~  N=rN  ~  N=i^*a  - 

Also,  for  Cosserat  surfaces  Cv,  the  expression  for  mechanical  power  cor- 

N 

responding  to  (3.22)  is 

P  =  Na  •  v  +  Z  •  w  +  Z  M^a  •  w., 


The  general  development  for  Cosserat  surfaces  Cv  outlined  above  is 

K. 

contained  in  a  paper  by  Green  and  Naghdi  (1976a,  Sec.  2)  which  deals  with 
application  of  the  theory  to  fluid  sheets  and  to  propagation  of  water 
waves.  When  K=l,  the  results  m  subsection  3.3  reduce  to  those  of 
subsection  3.2  for  a  Cosserat  surface  C. 


(3. 29) 


(3.30] 


4. 


Elastic  shells 


Within  the  scope  of  the  theory  of  a  Cosserat  surface  C  outlined  in 
Sec.  3,  we  discuss  briefly  the  constitutive  equations  for  elastic  shells 
in  the  presence  of  finite  deformation.  Preliminary  to  the  discussion  that 
follows,  we  assume  the  existence  of  a  strain  energy  or  stored  energy  per 

CL 

unit  mass  ^  =  ,t)  such  that  pip  is  equal  to  the  mechanical  power  defined 

by  (3.22),  i.e.. 


P  =  pip 


(4.1  ) 


In  the  development  of  nonlinear  constitutive  equations  for  elastic 

shells,  we  assume  that  the  strain  energy  density  ip  at  each  material  point 

of  C  and  for  all  t  is  specified  by  a  response  function  which  depends  on 

ot 

r,d  and  their  partial  derivatives  with  respect  to  0  .  But  since  the  response 
function  must  remain  unaltered  under  superposed  rigid  body  translational 
displacement,  the  dependence  on  r  must  be  excluded.  Thus,  the  constitutive 
assumption  for  the  strain  energy  density  can  be  written  as 


=  r  (r  ,d,d  :  X)  (■ 

and  we  also  make  similar  constitutive  assumptions  for  In  these 

constitutive  equations,  which  represent  the  mechanical  response  of  the 
medium,  the  dependence  of  the  response  functions  on  the  local  geometrical 
properties  of  a  reference  state  and  material  inhomogeneity  is  indicated 
through  the  argument  X. 

A  general  development  of  various  aspects  of  constitutive  theory  of 
elastic  shells  based  on  assumptions  of  the  type  (4.2)  or  variants  thereof 
is  given  in  Naghdi  (1972,  Sec.  13).  In  the  rest  of  this  section,  we  limit  the 
discussion  to  an  elastic  shell  which  is  homogeneous  in  its  reference  con¬ 


figuration  and  suppose  also  that  the  dependence  of  the  response  functions 
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on  the  properties  of  the  reference  state  occurs  through  the  values  of  the 
kinematical  variables  in  the  reference  state  (Carroll  and  Naghdi  1972). 
Then,  in  place  of  (4.2),  we  have 


(4.3) 


with  similar  assumptions  for  After  substituting  (4.3)  into  (4.1), 

by  usual  techniques  we  obtain  the  following  forms  for  the  constitutive 
equations : 


„a  dip  ,  3^  34> 

5  ’  “yT-  •  S  =  p3d  •  f  "  p  M 

.,a  ~  ~  ,a 


along  with  the  restriction 


(4.4) 


r 

-  .a 


~,a 


d  x  M  +  d 
3d  ~  ,a 


-,a 


(4.5) 


which  is  obtained  from  the  conservation  of  moment  of  momentum  and  which 
must  be  satisfied  by  the  response  function  ip  (Naghdi  1972,  Sec.  8). 

We  do  not  discuss  here  the  reduced  forms  of  the  above  constitutive 
equations  resulting  from  invariance  requirements  under  superposed  rigid 
body  motions,  but  for  such  reductions  refer  the  reader  to  Naghdi  (1972, 

Sec.  13).  Just  as  with  the  equations  of  motion,  it  is  necessary  in 
applications  to  specific  problems  to  obtain  alternative  forms  of  the  above 
constitutive  equations  or  their  reduced  forms  in  terms  of  tensor  components. 
Such  component  forms  may  be  expressed  with  respect  to  bases  a^,  or  d^,  or 
corresponding  bases  in  the  reference  configuration.  Reduced  forms  of 
(4.4)  have  been  utilized  extensively  in  Chapters  D  and  E  of  Naghdi  (1972). 
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5. 


Identification  of  the  assigned  fields  and  the  inertia  coefficients 


I 

! 


The  local  field  equations  (3.20)  in  the  mechanical  theory  of  a  Cosserat 
surface  have  the  same  forms  as  those  that  can  be  derived  from  the  three- 
dimensional  field  equations  (2.9)^  23^  suitable  integration  between  the 
limits  [recall  (2.17)  and  the  definition  of  a  shell-like  body  in 

section  2]  and  in  terms  of  certain  definitions  for  integrated  mass  density 
and  resultants  of  stress  (for  details,  see  Naghdi  1972,  Secs.  11-12 
or  Naghdi  (1974).  Similarly,  the  energy  equation  (3._’lj  has  the  same  form 
as  the  one  that  can  be  derived  from  the  energy  equation  in  the  three- 
dimensional  theory  by  suitable  integration  between  the  limits  £  anc^ 
in  terms  of  certain  definitions  for  integrated  internal  energy  density 


(5.5) 


,  _  .2  11, _  .2  22  33  .  12  „  13  ^  23, A 

ft.l)  '  8  ^1 ,2  8  +g  +  2  ^1 , 1^1 ,28  "^1,18  ~^1 , 28  ^ 


which  involves  the  partial  derivatives  of  5  (0a)  and  the  components  of  the 
metric  tensor  in  (2.3).  The  expression  for  f ^  can  be  stated  analogously. 
If  we  now  adopt  the  approximation  (2.25),  then  there  is  a  1-1  cor¬ 


respondence  between  the  two-dimensional  field  equations  that  follow  from 
the  conservation  laws  or  a  Cosserat  surface  and  those  that  can  be  derived 


from  (2.9)  -  ,  provided  we  identify  r  and  d  in  (2.25)  with  (3.1).  and 

(3.1)2,  respectively,  and  adopt  the  definitions  (5.1)-(5.4),  as  well  as  the 

definitions  of  the  resultants  mentioned  above.  A  similar  1-1  correspondence 

can  be  shown  to  hold  between  the  two-dimensional  energy  equation  in  the 

theory  of  a  Cosserat  surface  and  an  integrated  energy  equation  derived  from 

the  three-dimensional  energy  equation. 

The  various  quantities  in  (3.12)  are  free  to  he  specified 

in  a  manner  which  depends  on  the  particular  application  in  mind  and,  in 

the  context  of  the  theory  of  a  Cosserat  surface,  the  inertia  coefficients 
1  2 

y  ,y  and  the  mass  density  p  require  constitutive  equations.  Indeed, 
f  ,  ant*  V  3S  wel*  as  *b’^b  anc*  rb’  Can  identified  with  corresponding 
expressions  in  a  derivation  from  the  three-dimensional  equations  (for 

1  "> 

details,  see  Naghdi  1972,1979a).  Likewise,  p  and  the  coefficients  y  ,y~ 
may  be  identified  with  easily  accessible  results  from  the  three-dimensional 


theory. 

In  what  follows,  we  assume  that  the  above  identifications  have  been 

1  2 

made  and  that  the  quantities  p,y  ,y“,f^,H^  are  known  or  specified.  The 
knowledge  of  fc>%r  depends  on  the  nature  of  the  boundary  conditions  on  the 
major  surfaces  of  the  particular  shell-like  body  under  consideration:  they 
may  be  specified  as  known  quantities  on  the  surfaces  (2.20) ^  0  or  they  are 
unknown  (possibly  on  one  of  the  two  surfaces  (2.20) j  0  only)  and  must  be 
determined  as  part  of  the  solution  of  the  problem. 
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() .  Constrained  theories  of  shells 

A  development  of  a  constrained  directed  medium  in  the  3-dimensional 
theory,  with  particular  reference  to  an  incompressible  liquid  crystal 
having  a  single  director  of  constant  length,  is  contained  in  a  paper  of 
Green,  Naghdi  and  Trapp  (1970,  Sec.  6).  For  a  Cosserat  surface  with  a 
single  director,  a  number  of  constrained  theories  have  been  discussed  pre¬ 
viously.  These  pertain  to  a  class  of  shell-like  bodies  for  which  the 
director  is  constrained  to  be  of  constant  length  (Green  and  Naghdi  1974), 
an  incompressible  Cosserat  surface  (Green.  Laws  and  Naghdi  1974,  Green  and 
Naghdi  1976a)  and  a  class  of  fluid  sheets  in  which  the  director  is  constrained 
to  remain  always  parallel  to  a  fixed  direction  (Green  and  Naghdi  1977). 

A  special  case  of  the  constrained  theory  of  elastic  shells  discussed 

by  Green  and  Naghdi  (1974)  includes  that  for  which  the  director  is 

★ 

coincident  with  the  unit  normal  a^  to  the  surface  4.  This  special  form 
of  the  theory  can  be  brought  into  1-1  correspondence  with  that  of  a 
restricted  theory  of  elastic  shells  given  by  Naghdi  (1972,  Secs.  10,  15), 
where  the  director  is  not  admitted  and  the  basic  kinematical  ingredients 
that  occur  in  the  argument  of  the  strain  energy  response  function  are  a^ 
and  a^  ^  (compare  with  (4.3)).  Related  developments  include  the  construc¬ 
tion  of  a  theory  of  a  deformable  surface  with  simple  force  multipoles  by 
Balaban,  GTeen  and  Naghdi  (1967),  where  the  position  vector  r  and  its 

first  and  second  gradient  (r  ,r  .)  are  taken  as  the  basic  kinematic 

- ,  a  ~ ,  ap 

variables.  A  similar  theory,  but  less  general  than  that  of  Balaban  et  id. 
(1967),  is  given  by  Cohen  and  DeSilva  (1966a, 1968) .  For  additional 
related  comments  see  Naghdi  (1972,  Sec.  10). 

★ 

The  equations  resulting  from  such  a  constrained  theory  of  elastic  shells 
in  which  d =  a^  correspond  to  those  which  can  be  obtained  from  a  deriva¬ 
tion  of  shell  theory  under  the  so-called  Kirchhoff-Love  assumption  (sec 
Naghdi  and  Nordgrcn  1963). 
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In  this  section,  we  begin  by  considering  a  class  of  constraints  which 


are  linear  relations  between  the  kinematic  variables 


y,a  ’  ~N  ’  *N,ct  (N  =  1 ,2,  '  ‘  ‘  ,K) 


(<> 


in  the  form  (Green  and  Naghdi  1976a) 

.  v  ♦  Z  B^-w  ♦  Z  c^b  *  ww  =°  (M  =  o,2, . . .  ,Q)  ,  (<> 

~,a  N=  1  ~  ~N  N=1  ~  ~  ,a 

where  are  vector  functions  of  a^,d^,d^  ^  only  and  do  not  depend 

explicitly  on  the  variables  (6.1).  We  assume  that  each  of  the  functions 
rv  ^  Net  * 

N^.k  ,M  are  determined  to  within  an  additive  constraint  response  so  that 


,  kN-iN*kN  ,  .  vl* +W1 


(6 


where  NP‘,kN.^“  are  specified  by  constitutive  equations  and 


N°  ,  k*  ,  , 


(6 


which  represent  the  response  due  to  constraints  are  arbitrary  functions  of 
0a,t  and  are  workless.  Thus,  recalling  the  expression  (3.30)  for  mechanical 
power,  we  set 


-rfl 

N  •  V  + 


^  _n 
Z  k” 

N=1 


~N  + 


M^a  •  w 


N=1 


N,a 


=  0 


(6 


for  all  values  of  the  variables  (6.1)  subject 
(6.2).  It  then  follows  that 


N01 


Q 

z 

M=o 


,Ma 

A  Pk 


Q 

z 

M=o 


rMN 
?  Px 


to  the  constraint  conditions 


q 

Z 

M=o 


„MNa 

C  Pk 


(« 


The  development  between  (6.2)- (6.6)  is  similar  to  that  for  mechanical 
constraints  in  the  3-dimensiona*  theory  (see  section  30  of  Truesdell  and 
Noll  l!K>r.).  For  a  corresponding  thermodynamical  theory  of  a  continuum 
in  the  presence  of  thermo-mechanical  constraints  see  Green,  Naghdi  and 
Trapp  (1970)  and  Green  and  Naghdi  (1977). 
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a 

where  PM=P^(®  ,  t),  (M  =  o,  1 , . . .  ,Q)  are  arbitrary  functions  which  play  the 
role  of  Lagrange  multipliers. 

In  the  rest  of  this  section,  we  illustrate  the  nature  of  constrained 
theories  with  reference  to  two  particular  kinematical  constraints.  One  of 
these  constraints  is  that  appropriate  for  incompressible  media  and  the  other 
pertains  to  a  restriction  on  the  director  in  the  theory  of  a  Cosserat 
surface  C  with  a  single  director. 


6 . 1  Incompressible  Cosserat  surfaces 

The  conditions  representing  approximately  the  (3-dimensional) 
incompressibility  condition  (2.10)  may  be  derived  with  the  use  of  the 
approximation  (2.22)j.  However,  in  the  interest  of  brevity,  we  confine 
attention  to  a  special  case  of  (2.22)^  for  N=1  given  by  (2.25).  Under 
the  approximation  (2.25),  the  base  vectors  are  given  by  gQ  = aa  + 5d  a> 
g3  =  d,  where  a^  are  the  base  vectors  of  the  surface  5  =  0  calculated  from 
(2.25).  Then,  the  incompressibility  condition  (2.10)  may  be  expressed 
approximately  in  the  form 


3d  3d  3d  3d 


or  equivalently  as 


1 (d  •  a3)aa  -  (d  •  aa)a3  +  5(eaSd  g  *  d)]  •  v^a 

•  I»5  *  *!,«*  '•  *“*!,«  *  Ts1  •  ?  ♦  lie%  -  i  *  '  J.p 1  •  *.« - 

Ot 

where  in  (0.7)  and  (6.8)  use  is  made  of  the  notation  (3.29)  ,  and  e  ,  t: 

i ,  z  up 

denote  the  components  of  the  e-system  in  2-space  defined  by 

ap  -’-saB  11  22  1  2  21 

r.  =  a  e  ,  e  =  e  =  0  ,  e  =  -  c  =  1 


fcaB  '  a'emP  ’  C11  =  e22  °  ’  61 2  e21  ' 


We  now  generate  two  conditions  representing  incompressibility:  One  of  thes 


is  obtained  from  integration  of  (6.8)  with  respect  to  £  between  the  limits 


£j,£2  and  another  by  first  multiplying  (6.8)  by  £,  neglecting  terms 
3 

involving  £  and  then  integrating  the  resulting  equation  with  respect  to 
£  between  the  limits  £1(£9.  The  resulting  two  conditions  are: 

{y° [ (d  •  a  )aa  -  (d  •  aa)a  ]  +  y^^fd  *  d) }  •  v 

~  ~  ~  o  ~,P  -  ~,Ct 

+  {y°a3  +  £Cte[y1(ac(  x  d  >g)  +  %  y2  (d  >a  *  d  p) }  •  w 
+  {ea^  [y1  (ag  x d)  +  y2 (d  ^  * d) ] )  •  w  ^  =  0  , 

(y1  [ (d  *  a3)aa  -  (d  •  a“)a3]  +  yV^id  g  *  d)}  •  v  a 

+  (y1a3  + y2eae(aa  * d  g)}  ♦  w  + y2eaB(ag  *  d)  •  w  a  =  0  , 

K 

where  the  coefficients  y  are  defined  by 


(6.10) 


yK  =  /’2  £Kd£  ,  (K  =  o,l ,2) 


(6.11 ) 


It  is  perhaps  interesting  to  observe  that  in  special  circumstances  in  which 

* 

the  quantity  A  in  (5.1)- (5.2)  is  or  can  be  approximated  to  be  independent 

1  2 

of  £ ,  then  the  coefficients  y  and  y  in  (6.10)  will  have  the  same  numerical 

1  7 

values  as  the  inertia  coefficients  y  and  y~ ,  respectively. 

For  an  incompressible  Cosserat  surface  under  discussion,  from  (6.2) 

the  constraint  conditions  are 

.oa  t  Dol  A  „ola  .. 

A  •  v  +  B  •  w  +  C  •  w  =0 

~,a  ~  -  -  ~,a 

.la  ^  DH  ,  ~1  la 

A  •  v  +B  •  w  +  C  •  w  =0  , 

-  ,a  ~  ~  -  -  ,a 

and  the  corresponding  constrained  response  obtained  from  (6.6)  has  the  form 


(6.12) 


f  =-  (PoA^  +  PiA101)  , 

h  = '  (p0fa  +  Pi5n)  > 

f  *-(Pofla*P1clla)  . 


(6.15) 


where  F0-Pj  arc  the  Lagrange  multipliers.  Guided  by  the  two  conditions 

which  follow  from  (6.10),  we  select  the  vectors  A°a ,B°^ .C0*01 . which 

occur  in  (6.12)  to  have  the  special  values 


AKa  =  coeff.  of  v  in  (6.10) 

—  ~ 

v  1 

B  =  coeff.  of  w  in  (6.10) 

K  1  n 

C  =  coeff.  of  w  in  (6.10) 
~*a 


(K  =  o.l) 


Then,  it  follows  from  (b.1.3)  and  (b.14)  that  the  expressions  for  the  con- 

* 

straint  response  are  given  by 


('().  14  ) 


-rfi  „  (  ,  .  .a  ,  ,  a.  ,  ,,  a3  ,  „  . 

N  =-PQ[(d*a3)a  -  (d  •  a  )  a^j  -  dg*d  , 


k  =  -  P  a,  -  P,KaBa  xd  , 

o~3  1  ...  ,P 


'■z  Y*"p  f-a^d  x  d  , 

‘o  -  ,u  ...6 

2  fa , 


(o.l 3 ) 


rju  ,,  ah  .  z  pa  . 

M  =  -  P  e  a  x  d  -  y  p  e  d  x  d  . 

1  ~b  -  ‘o  ~ £ 

Che  arbitrary  coefficient  functions  Pq,Pj  are  related  to  the  Lagrange  multipliers 

•> 

i>  ,p,  and  v“p  can  he  expressed  in  terms  of  P  ,P,  as  follows: 

•  o  1  1  1  o  1  o  1 


=  Y  p  +  i  p.  ,  P. 
o  1  o  1  1  1  1 


1  2 

‘I  l\  +  Y  P, 


2  1 
Y  P  -Y  P, 


o  o2  .  1,2  ' 

Y  Y  -  ( Y  ) 


(  (,.  lb) 


In  obtaining  the  results  (6.10)  and  (b.lS),  no  identification  has  been 
made  between  the  surfaces  (3.1)^  in  the  theory  of  a  Cosserat  surface  C 
and  an  appropriate  reference  surface  in  the  (3-dimensional)  shell-like  body. 
Indeed,  different  values  for  the  coefficients  y^  in  (6.10)  will 
result  depending  on  the  choice  of  the  identification  with  the  surface 
f.  =  0.  For  example,  if  this  surface  is  chosen  between  the  major  surfaces 
of  the  shell-like  body  in  such  a  way  that  fj  =  then  the  coef¬ 

ficients  y  in  (5.S)  and  P  ,P,  in  (6.15)  become 

o  1 


o 

Y 


1  , 


1 


Y 


0 


2 


Y 


12 


P  = 
o 


(6.17) 


* 

Although  the  expressions  (6.15)  have  the  same  form  as  those  given  previously 
(Oreen  and  Naghdi  1976a,  hqs.  (4.2)),  they  are  not  the  same  in  view  of  the 
re  la  t  i ons  ( (> .  1  <> ) 
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and  the  incompressibility  conditions  (6.10)  reduce  to  those  used  by  Green 
and  Naghdi  (1976a,  Eqs.  (4.3))  for  a  directed  fluid  sheet  with  a  single 


director.  On  the  other  hand,  if  we  identify  (3.1)^  with  the  bottom  surface 

of  the  shell- like  body  so  that  £  =  0,  £0  =  1,  then  the  coefficients  y  and 

P  . P ,  become 
o  1 

o  ,  1  2  1 

Y=1,Y=%>Y=3» 

Po  =  P0  +  '*Pl  ’  Pl=^o  +  lPl  • 

For  a  complete  theory  of  an  incompressible  Cosserat  surface,  const itu- 

^  Ct  A 

tive  equations  are  required  for  the  quantities  N  ,k  and  M  but  a  discussion 
of  these  can  be  carried  out  as  in  Sec.  4. 


6.2  A  constrained  theory  with  director  along  the  normal  to  the  surface  of  (' 
We  turn  now  to  the  development  of  a  constrained  theory  of  a  Cosserat 
surface  in  which  the  director  is  always  along  the  normal  to  the  material 
surface  so  that 


((,  .  IK) 


d  *  aa  =  0  ,  d  =  <j>a3  ,  <f>  =  <J>  (0 T ,  t) 


(6.19) 


Differentiating  the  constraint  condition  (6.19)^  with  respect  to  time  and 


using  (3.3)  and  (3.4)  we  obtain 


d  •  v  +  a  •  w  =  0  , 

~  ~,a  ~a  ~ 


(6.20) 


which  represents  two  constraint  conditions.  From  (6.19)  ,  along  with  the  use 
of  (2.13)j  and  (2.15)  ,  follows  the  expression 


which  is  symmetric  in  a, 8-  Hence 


*. 


((>.22) 


C  (d,a’^  =  0  * 

D 

where  t  is  defined  in  (6.9).  Differentiating  (6.22)  with  respect  to  t  and 
*ct6 

observing  that  c  '  (d  a  •  a^)=0  in  view  of  (6.22)  and  the  fact  that 

'ah  up  .  , 

=  a  e  ,  we  also  have 


£  (d  •  v  0  +  a0  •  Vi  )  =  0 

-  ,a  ~,B  ~B  - ,a 


(6.23) 


as  a  third  constraint  condition. 

The  two  conditions  (6.20).  1  can  be  regarded  as  a  special  case  ((>.2) 

I,*- 

for  K=  1  and  with  the  coefficient  of  w  equal  to  zero.  Similarly,  (6.23) 

ft1 

is  a  special  case  of  (6.2)  for  M  =  o,  K  =  1  and  with  the  coefficient  of  w  equal 
to  zero.  Thus,  each  of  the  three  constraint  conditions  (6.20) ^  9  and  (6.23) 
may  be  viewed  as  a  special  case  of  (6.2)  with  coefficient  functions 
A  ,B  ,C  ,A  ,  B  ,C  conveniently  identified  as 


,ou  Bet,  ,.o  „ola  aB 

A  =  £  d,b  >  ?  =  e  >  9  =  c  > 

.Mu  .c» .  ..Ml  rMla  ...  , 

A  -  i-  d  ,  B  =  aM  ,  C  =o  (M  =1,2) 


(6.24) 


Now  according  to  (6.6)  and  with  the  help  of  (6.24),  the  expressions  for  the 
constraint  response  are 


rtx  ,  o  3ot,  r  Mxct , ,  ,  o  Bot .  a,, 

N  =  -  [p  e  d  g  ♦  f  p  6 Md J  =  -  [p  g  +  p  d]  , 

•"v 

r  v  M  a  rja  o  af 


(6.29) 


where  p  ,p  are  the  Lagrangian  multipliers,  and  in  line  with  the  notation  of 
(3.29)  and  that  of  subsection  3.2  we  have  set  k*  =  k,  M^a  =  n/* . 

In  anticipation  of  the  final  form  of  the  equations  of  the  constrained 
theory,  we  could  set  the  skew-symmetric  parts  of  both  and  equal  to 
zero  and  thus  require  also  the  vanishing  of  the  skew-symmetric  part  of 
which  is  equivalent  to  setting  p° = n .  However,  we  postpone  such  stipulations 
until  later  in  this  section,  and  retain  in  (6.25)  the  Lagrange  multiplier  p° 
which  arises  from  the  constraint  condition  (6.23). 
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Before  recording  the  modified  equations  of  motion  appropriate  for  the 

QL  QC  Y 

constrained  theory  under  discussion,  we  introduce  the  functions  S  =  S  (0  ,t) 


defined  by 


„a  r  a  a?  o  , 
S  =-  tP  Pjg] 


(6.2o) 


and  note  that 


ca  - 

la  13  la 


(6.27) 


Also,  for  convenience,  we  introduce  the  abbreviation 


f  =  f  -  (v  +  y*w)  ,  Z  =  Z-  (y*v  +  y2w) 


(6.28) 


Then,  after  substituting  (6.3)  and  (6.25)  into  (3.20)  and  making  use  of 

O  ,  ‘I  ,  J 


(6.26)  and  (6.27)  we  obtain 


N“,-pf  =-(*s“;3)|a  , 


ft*  - k-  *  0J  "  S\  ■ 


(6.20) 


a  x  Na  +  d  x  k  ♦  d  x  Ma  =  0  , 

-a  ~  ~  ~  ~,a  ~  ’ 

as  the  equations  of  motion  of  the  constrained  theory.  It  should  be  noted 

that  the  above  equations  involve  only  two  arbitrary  functions  of  position 

and  time  related  to  the  three  multipliers  p°,pa  by  (6.26).  Moreover,  (6.29) 

ot 

and  the  normal  component  of  (6.29)2  are  free  from  S  . 

A  further  reduction  of  the  system  of  equations  (6.29)  may  be  effected 
ot 

by  eliminating  S  between  (6.29) ^  ^ .  For  this  purpose  it  is  convenient  to 

refer  the  various  vector  quantities  in  (6.29)  to  the  base  vectors  a.  and 

★ 

write  the  equations  of  motion  in  tensor  components.  Thus,  we  write 


As  noted  earlier,  the  order  of  indices  in  (6.30)j  3  are  opposite  to  those 
used  in  Naghdi  (1972)  and  most  of  the  earlier  papers  on  the  subject. 


38. 


with  similar  expressions  for  hr.k.Nr.f ,H.  Recalling  (6 . 1 9) x  and  making  use 

o 

of  formulas  of  the  type  (6.30),  from  the  scalar  product  of  (6.29)^  with  aP 


and  again  with  a  we  deduce 


w^-*b?^Y)  =  0  ’ 

N3a-<(,ka-4)baM3Y-<|.  MaY  =  0  , 

Y  ,Y 


where  e  is  defined  in  (6.9). 

It  is  instructive  at  this  point  to  express  the  mechanical  power  (3.22) 
in  terms  of  the  tensor  components  (6.30).  To  this  end,  we  first  note  from 
(6.19)  and  (2.13)  -  that  the  tensor  components  of  w  and  w  referred  to  a 

A  I-  ,  Z.  yQL  ~ 


w*aa  =  -  <|>(v>a  •  a3)  ,  w  •  a3  =  <0 

w  •  a„  =  -  (4>b0  )  +  4>b  (v  a  •  aY)  -  <p  (v  Q  •  a,)  , 

~ ,a  3a  ay  ~,p  -  ,a 

yg 

w  •  a  ,  =  <f>  -  4>b  „  (v  •  a,)a' 

~,a  ~3  r,a  Y  aB  ~,Y  -3  ~ 

Then,  remembering  (6.4),  (6.5)  for  N=1  and  the  notation  (3.29).  7,  we  may 
write  (3.22)  as 

P  =  (N^a  -  <|>b^Y)  (v  a  •  ag)  +  k 5i  -  ^“(^b^)  ♦  M3a<j);a 
♦  -a3)  . 


Since  the  coefficient  of  (v  a  *  a3)  vanishes  identically  in  view  of 
(6.32)7,  the  last  expression  reduces  to 


(6.34) 


P  -  (N010  -  ^M“Y)  (v  a  •  ag)  +  k34> 

-  M^^FT)  ♦  „ 


~3a  ''a 

and  does  not  involve  the  components  N  and  k  .  Next,  with  the  help  of 

4>ka  +  4>M3Yba  =  N301  -  nPY<)>  which  follows  from  ( 6.32 )„,  the  component  form 
Y  »Y  *- 

of  the  equations  of  motion  (6.29) ^  ^  referred  to  a1  can  be  written  as 

*£a  .  gx,3a  3  n  M3a  ,  "3  n 

N  |  a '  baN  +P^  =  0  *  N  |ct  +  bfktN  +  pf  =  0  ’ 

(*M0a)|o  +  P^  =  n3°  ’  M3a|a  +  bBa^°‘'''3  +  P^  =  °  ’ 


(6.35) 


(6.36) 


where  in  recording  (6.35)  and  (6.36)  we  have  also  substituted  N  for  the 

quantity  (N'**  +  4>Sa) .  By  substitution  from  (6.36)^,  we  can  now  eliminate 

N301  from  (6.35)  .  In  this  way,  the  resulting  two  equations  may  be  put 

I  ,  z 

in  the  form 

{(if  4  ^  | a  *  ba(^Y  *  »“>  |Y}^  +  *  ?“)  |$a  +ba3(?  *  ^)}?3 

+  (pf  +  (p<j)iiaa3)  |a)  =0  .  (6.37) 

In  a  general  theory  of  an  elastic  Cosserat  surface  (Sec.  4),  constitutive 
equations  for  both  the  symmetric  and  the  skew- symmetric  parts  of  Na®,M°^  can 
be  provided  through  the  expression  for  mechanical  power.  Here,  however,  since 
b^g  is  symmetric,  the  term  -M^^b^)  in  (6.34)  provides  constitutive  equations 
for  only  the  symmetric  part  of  M01® .  Moreover,  the  quantity  (Na^  -  4>bbNPY)  is 
symmetric  by  virtue  of  (6.32) j  and  the  two  differential  equations  resulting  from 
(6.37)  involve  only  the  symmetric  parts  of  N0^  and  .  Thus,  in  line  with 
classical  results  in  shell  theory,  in  order  to  obtain  a  determinate  theory  we 


now  put 


,  =  %(M“6  -  M06) 


(6.38) 


In  summary,  the  relevant  system  of  equations  of  the  constrained  theory  under 

discussion  are  given  by  (6.37),  the  normal  component  of  (6.  29)2,  i-e.. 

ct6 

(6.36)2  and  the  skew-symmetric  part  N  is  determined  from  (6.32) j.  This 


Instead  of  introducing  (6.38)j,  in  anticipation  of  the  fact  that  makes 

no  contribution  to  the  mechanical  power,  at  the  outset  we  could  have  absorbed 
M‘ab*  into  M I J  or  equivalently  into  FP  in  (6.25). 


completes  the  development  of  the  constrained  theory  in  which  the  director 


is  constrained  to  have  the  form  (6.19)^. 

If  in  addition  to  (6.37)  we  also  set  the  multiplier  p° = 0,  then 
^[a8]  _q  ancj  hence  skew-symmetric  =0.  It  then  follows  that 


ca  a 
S  =  -  p 


N**  =  N**  -  *Sa  ,  ka  =  ka  -  SQ 


and  the  relevant  equations  of  motion  of  the  determinate  constrained  theory 
remain  as  before.  It  is  of  interest  to  examine  the  reduction  of  the  fore¬ 
going  development  when  <p  =  1 .  In  this  case,  we  have  d  =  instead  of 
(6.19)^  and  the  resulting  equations  are  identical  with  those  of  a  restricted 
theory  discussed  by  Naghdi  (1972,  Secs.  10  and  IS).  The  results  with  <}>  =  1 
can  also  be  brought  into  correspondence  with  a  special  case  of  the  con¬ 
strained  theory  discussed  by  Green  and  Naghdi  (1974)  or  those  contained  in 
the  paper  of  Naghdi  and  Nordgren  (1963) . 

The  nature  of  the  boundary  conditions  in  the  theory  of  a  Cosserat 
surface  C  discussed  in  subsection  3.2  is  clear  from  the  expression  for  the 

rate  of  work  R  of  contact  force  and  contact  director  force  over  the  closed 
c 

boundary  curve  3£P,  namely 

R  (£?)  =  J  (n  •  v  +  m  •  w)ds 
c  3ff~  ~  ~ 

However,  in  a  constrained  theory  of  the  type  discussed  in  subsection  0.2, 

the  question  of  the  boundary  conditions  must  be  reconsidered  in  view  of 

* 

the  reduction  in  the  number  of  differential  equations  .  Since  the  develop¬ 
ment  of  the  reduced  boundary  conditions  is  similar  to  that  of  a  restricted 

The  number  of  the  relevant  scalar  differential  equations  of  the  constrained 
theory  is  five  as  compared  to  the  nine  scalar  equations  in  the  theory  of 
subsection  3.2. 


theory  (Naghdi  1972,  p.  552),  our  discussion  will  be  brief. 
Recalling  (6.30)^  ^  3  an^  (6.33),  from  (6.41)  we  obtain 

R  (29  =  /  v  { (NyCX  -  4>b^MSa)v  +  N^v  +  M^i  -  <j>MyV  }ds 


(6.42) 


Let  3/3v  stand  for  the  directional  derivative  along  the  unit  normal  v  to 
the  boundary  curve  3£?  and  let  3/3s  denote  the  directional  derivative  along 
the  tangent  to  3 £P.  Then,  provided  the  quantities  in  (6.42)  are  single¬ 
valued  on  a  (sufficiently  smooth)  closed  curve  3£P,  with  the  use  of 


8v3  8v3  6 

v3,y  =  IkT  VY  '  IT  £ybv 


(6.43) 


and  an  integration  by  parts,  (6.42)  can  be  reduced  to 


-  UP VpV3’G  fcT+H*}d8  > 


(6.44) 


where 


p6  =  (N001  -  4>beMYa)v  ,  G  =  4>My“v  v  =  d>M('y0t') v  v 
f  y  a  y  a  y  a 

3  ...3a  3  ,,,_y  a  _  ,,  ,,3a 

p  =  4»N  v  -  (4>M '  v  e  „vM)  ,  H  =  M  v 

v  a  3s  a  yp  a 


(6.45) 


The  nature  of  the  reduced  boundary  conditions  of  the  constrained  theory  is 
now  clear  from  (6.44)  and  (6.45). 
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7.  Additional  remarks  on  shells 

The  theory  of  Cosserat  surfaces  can  easily  allow  for  the  effect  of 
surface  tension  (Naghdi  1972,  p.  547  and  1974)  and  can  accommodate  the 
specification  of  either  tractions  or  displacements  on  major  surfaces  of  the 
shell-like  (or  sheet-like)  bodies  for  application  to  various  interfacial  and 
contact  problems.  Even  the  theory  of  a  Cosserat  surface  with  a  single 
director  can  be  used  to  formulate  a  fairly  broad  class  of  contact  problems 
of  elastic  shells  and  plates,  as  discussed  by  Naghdi  (1975a).  The  relevance 
and  applicability  of  the  basic  theory  of  a  Cosserat  surface  to  problems  of 
an  incompressible,  inviscid  fluid  sheet  is  discussed  by  Green,  Laws  and 
Naghdi  (1974)  and  by  Green  and  Naghdi  (1975,1976a).  The  nonlinear  dif¬ 
ferential  equations  derived  in  these  papers  include  the  effects  of  gravity 
and  surface  tension  and  are  also  valid  for  propagation  of  fairly  long  water 
waves  in  a  stream  of  initial  variable  depth.  A  discussion  of  an  incompres¬ 
sible  viscous  fluid  sheet,  along  with  further  recent  developments  on  the 
subject,  can  be  found  in  the  papers  of  Green  and  Naghdi  (1976a, b ; 1977a ; 1979c) . 
The  basic  theory  is  also  applicable  to  problems  of  cell  membranes,  as  has 
been  emphasized  by  Ericksen  (1979). 

In  the  remainder  of  this  section  we  briefly  comment  on  some  special  cases 
of  the  general  theory  and  also  mention  some  recent  researches  which  bear  on 
the  various  aspects  of  elastic  shells.  Although  these  developments  will  be 
described  mainly  in  the  context  of  a  mechanical  theory,  some  recent  results 
pertaining  to  thermal  effects  in  shells  are  also  discussed. 

The  well-known  membrane  theory  of  shells  can  be  obtained  as  a  special 
case  of  the  general  theory  by  essentially  suppressing  the  effect  of  the 
director  and  corresponding  kinetical  variables  and  this  is  discussed  briefly 
in  Naghdi  (1972,  Sec.  14).  A  development  of  another  special  theory,  known 
as  the  inextensionai  theory,  wherein  the  length  of  each  clement  of  the 
surface  of  4  is  assumed  to  remain  constant  throughout  all  motions  is  also 


» 


contained  in  Naghdi  (1972,  Sec.  14).  Similarly,  a  nonlinear  restricted 
theory  of  shells  by  direct  approach,  motivated  mainly  by  the  classical 
theory  corresponding  to  Kirchhoff-Love  theory  of  shells  and  plates,  is 
given  by  Naghdi  (1972,  Secs.  10  and  15).  Related  constrained  theories  of 
an  elastic  Cosserat  surface  are  already  mentioned  in  Sec.  6  and  need  not 
be  repeated  here. 

The  nonlinear  constitutive  equations  in  Sec.  4  are  valid  for  an  elastic 
Cosserat  surface  which  may  be  anisotropic  with  reference  to  preferred  direc¬ 
tions  associated  with  material  points  of  <5.  A  general  discussion  of  material 
symmetries  for  shells  is  given  by  Naghdi  (1972,  Sec.  13).  Carroll  and  Naghdi 
(1972)  have  subsequently  examined  the  influence  of  the  reference  geometry  on 
the  response  of  elastic  shells  by  assuming  the  existence  of  a  local  preferred 
state  of  the  body  and  then  stipulating  that  the  influence  of  the  reference 
geometry,  as  in  (4.3),  occurs  through  the  values  of  the  constitutive  variables 
in  the  preferred  state.  Material  symmetry  restrictions  for  elastic  shells 
have  been  discussed  also,  from  a  different  point  of  view,  by  Ericksen  (1972a, 
1973b)  who  has  also  indicated  (Ericksen  1973)  a  comparison  with  the  results 
contained  in  the  paper  of  Carroll  and  Naghdi  (1972). 

Some  general  aspects  of  wave  propagation  in  elastic  shells,  based  on 
the  theory  of  a  Cosserat  surface  have  been  discussed  by  Ericksen  (1971).  A 
related  study  on  the  subject,  limited  only  to  wave  propagation  in  a  surface 
not  endowed  with  a  director,  was  given  earlier  by  Cohen  and  Suh  (1970).  The 
theory  of  small  deformation  superposed  on  a  large  deformation  of  an  elastic 
Cosserat  surface,  along  with  related  problems  of  stability  and  vibrations  of 
initially  deformed  plates,  is  discussed  by  Green  and  Naghdi  (1971).  Related 
developments  concerning  plane  waves  and  stability  of  elastic  plates  are  given 
by  Ericksen  ( 1973c , 1974) .  For  a  system  of  linear  equations  characterizing  the 
initial  mixed  boundary-value  problem  of  elastic  shells,  Naghdi  and  Trapp  (1972) 
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have  obtained  a  uniqueness  theorem  without  the  use  of  definiteness  assumption 
for  the  strain  energy  density.  This  result  (Naghdi  and  Trapp  1972)  holds 
for  nonhomogeneous  and  anisotropic  shells  undergoing  small  motions  super¬ 
posed  on  a  large  deformation. 

In  still  another  study,  the  theory  of  a  Cosserat  surface  has  been 
employed  by  Naghdi  (1975a)  to  formulate  contact  problems  of  shells  and 
plates  mentioned  above.  In  the  derivation  of  shell  theory  from  the  3- 
dimensional  equations,  equations  of  motion  in  terms  of  resultants  and 
detailed  consideration  of  constitutive  equations  for  shells  are  usually 
obtained  relative  to  an  interior  surface,  rather  than  one  of  the  major 
surfaces  of  the  shell-like  body  which  may  be  the  contacting  surface;  the 
interior  surface  ordinarily  is  identified  with  the  middle  surface  of  the 
shell  or  plate  in  the  reference  configuration.  In  the  development  of  shell 
theory  by  direct  approach,  although  the  material  surface  of  S  may  be 
identified  with  any  surface  of  the  (3-dimensional)  shell-like  body, 
nevertheless  the  complete  discussion  of  constitutive  equations  and  the 
identification  of  the  inertia  coefficients  and  the  assigned  fields  may 
again  require  explicit  use  of  a  reference  surface  in  the  shell-like  body. 

For  certain  problems  it  is  more  natural  and  conceptually  more  appealing  to 
select  one  of  the  two  major  surfaces  as  the  reference  surface  but  then  the 
detailed  available  development  of  the  constitutive  equations,  as  well  as 
identification  of  such  quantities  as  the  inertia  coefficients,  have  to  be 
reconsidered  relative  to  the  new  surface.  This  problem  can  be  resolved  by 
deriving  appropriate  transformation  relations  (Naghdi  1975a),  which  relate 
the  kinetic  variables  n,k,m  (and  hence  the  response  functions)  in  the  two 
formulations.  The  results  (Naghdi  1975a)  are  applicable  to  any  shell-like 
medium  and  their  validity  is  not  limited  to  elastic  shells  alone. 

Controllable  solutions  in  the  theory  of  a  Cosserat  surface  have  been 
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studied  by  Crochet  and  Naghdi  (1969),  Ericksen  (1972b)  and  Naghdi  (1975b). 

In  a  more  recent  study,  Naghdi  and  Tang  (1977)  have  discussed  controllable 
deformations  that  can  be  maintained,  in  the  absence  of  body  force,  in  every 
isotropic  elastic  membrane  by  the  application  of  edge  loads  and/or  uniform 
normal  surface  loads  on  the  major  surfaces  of  the  thin  shell-like  body. 

The  static  solutions  of  finitely  deformed  membranes,  which  are  valid  for 
both  compressible  and  incompressible  materials,  are  obtained  with  the  use 
of  a  strain  energy  response  function  which  depends  on  the  metric  tensor  of 
the  membrane  it  its  deformed  configuration.  The  main  results  are  summarized 
by  several  theorems  and  their  corollaries  in  accordance  with  three  mutually 
exclusive  cases  for  which  the  initial  undeformed  surface  of  the  membrane 
(which  may  be  a  sector  of  a  complete  or  closed  surface)  is,  respectively, 
developable,  spherical  and  a  surface  of  variable  Gaussian  curvature 
satisfying  certain  differential  criteria.  The  corresponding  deformed 
surfaces  are,  respectively,  a  plane  or  a  right  circular  cylinder,  a  sphere 
and  a  surface  of  constant  mean  curvature.  These  results  are  exhaustive  in 
that  they  represent  all  finite  deformation  solutions  possible  in  every 
isotropic  elastic  material  characterized  by  the  strain  energy  response 
mentioned  above.  Also  discussed  in  the  paper  of  Naghdi  and  Tang  (1977)  are 
some  special  cases  of  the  general  results  and  several  families  of  solutions 
in  terms  of  an  alternative  description  which  should  be  useful  in  application 
and  which  permit  easy  interpretation. 

The  development  of  the  theory  of  Cosserat  surfaces  in  Sec.  5  i - 
carried  out  within  the  scope  of  the  purely  mechanical  theor\ .  In  earlier 
work  on  thermo-mechanical  theory  of  shells  by  direct  approach  (i.reen  and 
Naghdi  1970,  Naghdi  1972),  only  one  temperature  field  was  admitted  and  tins 
allowed  for  the  characterization  of  temperature  changes  along  some  reference 
surface,  such  as  the  middle  surface,  of  the  (3-dimensional)  shell-1  ike  body. 


but  not  for  temperature  changes  along  the  shell  thickness.  The  latter 
effect  has  been  incorporated  recently  by  Green  and  Naghdi  (1979a)  into  the 
thermo-mechanical  theory  of  Cosserat  surfaces,  together  with  appropriate 
thermodynamical  restrictions  arising  from  the  second  law  of  thermodynamics 
for  shells. 


i 


I 
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8. 


The  basic  equations  in  direct  notation 


For  some  purposes  it  is  convenient  to  have  available  the  basic  equa¬ 
tions  for  a  Cosserat  surface  in  a  direct  (coordinate- free)  notation  and 
this  is  the  main  purpose  of  the  present  section.  As  will  be  evident 
presently,  the  forms  of  the  basic  equations  in  coordinate- free  notation 
are  very  similar  to  those  of  the  corresponding  equations  in  the  classical 
3-dimensional  theory  and  thus  may  be  more  suitable  in  the  discussion  of 
general  theorems  or  in  developments  which  parallel  those  in  the  3-diinension.il 
theory. 

As  in  the  papers  of  Carroll  and  Naghdi  (1972)  and  Naghdi  (1977),  we 
introduce  the  notations  grad  and  Grad  to  denote  the  right  spatial  and 
material  gradient  operators,  respectively,  with  respect  to  the  position  on 
the  surface  4  in  the  current  configuration  and  on  the  surface  in  the 
reference  configuration.  The  corresponding  divergence  operators  will  be 

denoted  by  div  and  Div,  respectively.  In  particular,  for  a  vector-val ued 

a  * 

function  V(0  ,t),  we  write 


grad  V  =  Va®d“ 

Grad  V  =  V  ®Da 

-,a  ~ 


div  V 

Div  V 


(8-1) 


where  the  symbol  ®  denotes  the  tensor  product.  Also,  the  spatial  surface 
gradient  and  the  spatial  surface  divergence  operators  are  defined  by 


grad  V  =  V  aa  ,  div  V  -  V  *a 
6  s  ,a~  s  ~  ~,ct  ~ 


(8.2) 


We  take  this  opportunity  to  correct  an  error  in  a  previous  paper  (Naghdi 
1977).  The  definitions  (2.9)j  2  of  Naghdi  (1977)  should  be  replaced  with 
those  in  (8.1)2  2  of  the  present  paper  with  da  defined  through  (3.5). 

Also,  the  "div"’operator  in  (3.10)2  of  Naghdi  (1977)  should  be  replaced  by 
”di vs”  in  (8 . 2) 2 .  The  definitions  ( 2 . 9) 3  4  of  Naghdi  (1977)  remain  unchanged 
since  previously  (Naghdi  1977)  the  director  in  the  reference  configuration 
was  specified  to  have  the  form  D=DA3.  Except  for  the  modifications  noted, 
all  other  results  in  the  paper  of  Naghdi  (1977)  remain  intact. 


for  all  scalar-valued  functions  V  and  all  vector-valued  functions  V. 

We  introduce  a  measure  of  deformation  by  the  tensor  F,  namely^ 

F  =  d.  O  D1  =  Grad  r + d7  0  D3  ,  (8.3) 

~i  ~  ~3  ~ 

and  in  view  of  the  notations  (3.5)  and  (3.9)  we  observe  that 


F  D 
-  ~a 


F  A. 


F  D3  =  F  D  =  d  =  dj 


(8.4) 


From  the  definition  of  the  determinant  of  a  second  order  tensor  T  given  by 


det  I^Yi^2~3^  =  II  ~1’I  ~2’~ 

for  all  arbitrary  vectors  vi>v2,v3>  the  conditions  (3.1)^  and  (3.7)^ 
we  obtain 


det  F  =  [d1d2d3]/[D1D2p3]  >  0 


(8.5) 


The  tensor  F,  a  linear  operator  on  vectors  in  3-space,  is  nonsingular;  and 
there  exists,  therefore,  the  inverse  deformation  gradient  tensor  F  1 
defined  by 

F"1  =  D.  0  d1  .  (8.0) 

~i 

The  inverse  operator  F  1  transforms  vectors  in  the  present  configuration  into 
vectors  in  the  reference  configuration,  i.e., 

F_1d.  =  D.  (8.7) 

-  ~i  ~i 


and  it  follows  that 


This  definition  of  F  is  the  same  as  that  used  by  Naghdi  (1977).  The 
symbol  F  in  the  paper  of  Carroll  and  Naghdi  (1972)  stands  for  a  different 
quantity.  The  term  Grad  r  in  (8.3)  corresponds  to  the  deformation  gradient 
tensor  F  in  the  paper  of  Carroll  and  Naghdi  (1972). 
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(8.8) 


F~ 1 F  =  F  F'1  =  I  =  d.Sd1  =  D.0D1  , 

~  ~  ~  ~  ~  ~i  -  ~i  ~ 

where  I  is  the  unit  tensor  in  3-space.  We  also  introduce  here  the  director 
gradient  tensor  G  by 


G  =  Grad  d  =  d,  ®d“  =  d  ®Da  . 

~3,Qi  ~  ~,a  ~ 

Recalling  the  definitions  (3.3)^  ^  for  the  velocity  and  the  director 

velocity  and  since  a  =  v  ,  we  have 
1  -a  ~,a* 

F  =  d.  ®  D1  =  d  ®Da  +  d,®  D3  =  v  ®Da  +  w®D3  , 

-l  ~  ~a  -  -3  ~  -,a  ~  ~ 

G  =  d,  ®Da  =  w  ®Da  . 

-  ~3,a  -  -,a  - 

Also, 


(8.10) 


F  F  1  =  di®di  =  grad  v  +  w®d3  , 


G  F  1  =  w  ®da  =  grad  w 

-  -  ~,a  ~  ~ 


(8.11) 


Having  disposed  of  the  main  kinematical  results  in  terms  of  the 
gradient  tensors  F,G  and  their  rates,  we  now  turn  to  kinetical 

quantities.  The  expressions  corresponding  to  ( 3 . 1 8 ) ^  for  the  contact 

5 

force  n  and  the  contact  director  m  can  now  be  expressed  in  the  form 


n  =  Nv  ,  m=Mv, 


(8.12) 


yThe  second  order  tensors  N,M  in  (8.12)  and  their  tensor  components  Nlu,Mlul  in 
(8.13)  are  the  trapsnose  of  the  corresnonding  quantities  in  Naghdi  (1977).  The 
components  Nia,  Mia  were  used  in  the  paper  of  Green,  Naghdi  and 
Wainwright  (1965)  but  subsequently  their  transpose,  namely  N011  and  Nf11  , 
were  adopted  in  subsequent  papers  so  that  the  notation  would  be  in  agree¬ 
ment  with  that  of  the  classical  shell  theory.  It  may  be  noted  that  in 
terms  of  the  latter  notation,  instead  of  (8.12),  one  would  have  n=N^y> 
m =  M^u,  where  the  superposed  T  denotes  transpose.  Compare  (3.6)  and  (3.10) 
of  Naghdi  (.1977)  with  (8.12)  and  (8.15)  of  the  present  paper. 
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with  the  second  order  tensors  N,M  defined  by 


i~)i  1  Ol  ( 1  ■  y 

N  =  N  *d  =  N  d.*d  ,  N  =  N  d  , 

-(x  -a  ~  ~  - 


. .  .  ,u ,  , ,  i  a  ,  „  ,  „  .a 

M  =  H  {(id  =  M  d .  W d  ,  M  =  M  d 

fit  -  l  ..a 


(  K .  1  S  ) 


ci  i 

which  also  relate  the  tensors  N,M  to  N  ,M  in  (4.8).  Also,  for  convenience, 
he  introduce  a  second  order  tensor  K  through 


K  =  k^d„  =  kJd.  «od„  ,  k  =  Kd3 

~  ~j>  ~i  ~3 


(8.14) 


With  the-  use  of  is. I  si  and  by  usual  procedures,  from  the  conservation  laws 
Id. I o  |  follow  the  local  equations 


p  +  p  di v  v  =  0  , 


di v  N  +  pf  =  p (v  +  y1 w)  , 


]  ‘  2' 

d  i  v  M  +  p  Z  -  k  =  p  ( v  v  +  v  w )  , 

s  -  -  ' 


(  8  .  I  S  ) 


[N  +  K  +  M(('i  I-'1)1]  =  [N  +  K  +  M(G  F-1)1]1 


which  are  equivalent  to  (A. 20).  Also,  by  the  definition  of  the  right 
divergence  of  a  tensor  field,  we  have 

div  N  =  Na.  ,  div  M  =  m\  .  (8.  !M 

s  ~  -  |  a  s  -  s  |  a 

It  is  interesting  that  the  last  statement  in  (8.  IS)  is  analogous  to  the  sxmiiietrx 
of  the  stress  tensor  in  the  S- d  imensi  ona  1  theory.  In  part  ivulai  ,  it  max  lie 

rt  it 

observed  that  a  *N  ,  d  <  k  and  d  *M  are,  respectively,  the  axial  vectors 

~'x  ,t 

of  |.\-N*|,  [  K  -  K 1  |  and  |M(C  l;  *)!  -M((l  1-  S|.  Furthermore,  in  terms  of  the 

Kinetical  quantities,  \,M,k  in  (  8  .  I  2  )  -  (  8  .  1  4  )  and  the  rate  quantities  ( ts  .  If)  j 
the  mechanical  power  becomes 

I’  t  r t  |  I-1  (N  +  K)  +  C i 1 M |  ( 1  1  l  1 }  18.1") 


With  reference  to  constitutive  equations  for  elastic  shell 
ot  the  Kineniatical  variables  in  section  1,  in  now  einplov  the  va 
and  (8.;)).  Thus,  corresponding  to  constitutive  assumption  M.S 
w  r  i  t  e 

if  -  '+>  ( !• .  (i  ;  |ji> )  , 

where 


,  i;  =  brad  li  -  0.  I' ' 
lb  -  a  .  .i 


alonit  with  similar  assumptions  for  \,k,M.  Then,  with  the  usi.  o 
and  (8.17),  by  usual  techniques  we  obtain  the  following  altcrn.i 
of  the  constitutive  equations: 


N  +  K 


the  first  of  which  can  be  resolved  into 


N  -  c  S  d>'*d  ) 

•I-  .  x 


.  K  --  ,  (D'^d_) 

■  a  I-  .  ..•> 


M  ai,  the  response  function  if  is  restricted  l>v 


S  -  S 


i  + t;1 


jm  i,  rji  8 „  ai,'1 

-  ( I)  d  )  +  -  r  ( H  *  d,)  +  b 


-  ,  i  list  cad 

1  ab  I  e-.  (  8  .  i  ! 

,  we  in  iw 

(  8  .  I  8  | 

IS. IS! 

'  (Mi 
i ve  forms 

(  8  .  J  ( 1 1 

(8.711 


S 


(8.77.1 


Part  B 


Plastic  rods:  A  direct  formulation 


In  Part  B  (Secs.  9-13),  we  first  summarize  the  main  kinematics  and 
tlie  basic  principles  of  the  theory  of  Cosserat  (or  directed)  curves  and  then 
discuss  the  constitutive  equations  for  elastic  rods,  as  well  as  some  related 
aspects  of  the  basic  theory  and  recent  developments  on  the  subject.  Although 
we  are  concerned  here  mainly  with  the  purely  mechanical  theory  involving 
appropriate  forms  of  the  conservation  laws  for  mass,  linear  momentum,  director 
momentum  and  moment  of  momentum,  we  also  include  a  statement  of  the  conserva¬ 
tion  of  energy.  The  latter  provides  motivation  in  the  development  of 
certain  constitutive  equations,  such  as  those  for  an  elastic  material,  and 
in  the  discussion  of  aspects  of  some  special  solutions  involving  jump  in 
energy.  The  contents  of  Part  B  are  as  follows: 

9.  The  basic  theory  of  Cosserat  curves 

9.1  Kinematics  of  a  Cosserat  curve  5?. 

9.2  Basic  principles  of  a  Cosserat  curve  31. 

9.3  Hierarchical  theories  of  Cosserat  curves. 

ill.  Plastic  rods 

11.  Identification  of  the  assigned  fields  and  the  inertia  coefficients 

12.  Additional  remarks  on  rods 

13.  Itie  basic  equations  for  clastic  rods  in  direct  notation 


y .  The  basic  theory  of  Cosserat  curves 

Having  defined  a  (three-dimens icmal j  rod-like  body  in  section  2,  we  now 

formally  introduce  a  direct  model  for  such  a  body.  Thus ,  deformable  media 

which  are  model  led  bv  a  material  curve  J'  embedded  in  a  fuel i dean  3- space, 

together  with  l.  (L>2)  deformable  vector  fields --  cal  led  directors  -- 

attached  to  every  point  of  the  material  curve  are  called  Cos sera t  curves  or 

directed  curves  and  may  he  conveniently  referred  to  as  SR  (k=  1,2,...,N). 

The  directors  which  are  not  necessarily  along  the  unit  principal  normal  and 

the  unit  binormal  vectors  to  the  curve  have,  in  particular,  the  property 

that  they  remain  unaltered  in  length  under  superposed  rigid  bode  motions. 

In  the  absence  of  the  directors,  we  merely  have  a  1 -dimens i ona 1  material 

curve  f  which  can  serve  as  a  model  for  the  construction  by  direct  approach 

of  string  theory.  The  relationship  between  the  number  of  directors  1. 

and  the  number  K  which  identifies  the  order  of  the  hierarchical  theory  of 

K 

Cosserat  curves  can  be  shown  to  be  L  =  E(N  +  1)  so  that  (see  Naghdi  19791') 

1 

L  =  K(K+3)/2  . 

With  K  =  1 ,  the  directed  curve  is  a  body  5?^  =51  comprising  a  material  curve 
anti  two  deformable  directors  attached  to  every  point  of  the  material  curve 
of  5?.  The  latter  is  the  simplest  model  for  the  construction  of  a  general 
bending  theory  of  slender  rods;  and,  for  simplicity,  we  restrict  attention 
to  this  particular  model  in  most  of  the  development  in  section  9. 

We  now  turn  to  a  brief  account  of  the  basic  theory  of  a  Cosserat  curve. 

<).  i  kinemat  ics  of  a  Cosserat  curve  5? 

I.et  the  particles  of  the  material  curve  1  of  fR  he  identified  bv  means 
of  the  convected  coordinate  l  and  let  the  curve  occupied  by  J'  in  the  present 

.1. 

\  brief  account  of  the  more  general  theory  for  Cosserat  curves  is  indicated 
at  the  end  of  this  section. 


configuration  of  ‘R  at  time  t  be  referred  to  as  £.  Let  r  and  d^  (u  =  1,2) 
denote  the  position,  vector  of  a  typical  point  of  £  and  the  directors  at  the 
same  point,  respectively,  and  also  designate  the  tangent  vector  to  the 
curve  £  by  a^.  Then,  a  motion  of  the  Cosserat  curve  is  defined  by  vector¬ 
valued  functions  which  assign  a  position  r  and  a  pair  of  directors  d  to 

★  + 

each  particle  of  at  each  instant  of  time,  i.e. 


5  =  r^’t}  ’  ~a  =  ’  ^1^2-a3*  >  ° 


where 


?3  =  =  at 


The  condition  (9.2)^  ensures  that  the  directors  d^  are  nowhere  tangent  to 
£  and  that  d^.d^,  never  change  their  relative  orientation  with  respect  to 
each  other  and  a,.  The  velocity  and  the  director  velocities  are  defined  by 


v  =  r  ,  w  =  d 
~  -a  ~a 


(9.1) 


and  from  (9.3)  and  (9.4)  we  have 


~3  ”  9£  ’ 

where  a  superposed  dot  denotes  material  time  differentiation  witli  respect  to 
t  holding  £  fixed. 

It  is  convenient  to  introduce  here  a  slightly  different  notation  than 
that  adopted  in  a  number  of  previous  papers,  e.g.,  Naghdi  (1979a).  Thus,  we 


lor  convenience,  we  adopt  the  notation  for  r  in  (2.1b)  and  (2.30)  also 
for  the  surface  (9.2)j.  This  permits  an  easy  identification  of  the  two 
curves,  if  desired.  The  choice  of  positive  sign  in  (9.2)^  is  for 
definiteness.  Alternatively,  it  will  suffice  to  assume  that  ldjd>a-.|/0 
with  the  understanding  that  in  any  given  motion  the  scalar  triple  product 


I  d  ,  d  ,a  ,  |  i  s  e  i  t  her  >0  or  •  0  . 
'  I  2  •> 


(9.0) 


^3  =  “3  ’  *i  =  (ix*?3> 


and  observe  that  in  view  of  (9.2)^  and  (9.6),  d^.d^d,  are  linearly 
independent  vectors.  Hence,  we  may  introduce  a  set  of  reciprocal  vectors 
d*  such  that 


d.  •  dJ  =  <5J  , 

-1  ~  1 


where  6.  is  the  Kronecker  symbol  in  3- space.  Whenever  desirable,  the 
notations  d.  =  (d  ,d^,d,)  and  (d  , a ,)  will  be  used  interchangeably  throughout 
Part  B  depending  on  the  particular  context.  Consider  now  a  reference  con¬ 
figuration,  not  necessarily  the  initial  configuration,  of  the  Cosserat 
curve  <H .  In  the  reference  configurat  ion,  let  the  material  curve  of  IK  be 
referred  to  by  X  and  designate  the  unit  principal  normal,  the  unit  binomial 
and  the  tangent  vector  to  X  by  A  ,A  and  A.,  respectively.  further,  let 
R  and  1)^  (a  -  1,2)  stand  for  the  position  of  a  typical  point  of  X.  anil  the 
directors  at  the  same  point,  respectively.  Then,  in  the  reference  con¬ 
figuration  we  have 


=  R(fJ  ’  !?a  =  >  0 


where 


*3  =  *3^  -  t 

and  (9.8)^  ensures  that  are  nowhere  tangent  to  the  curve  X^.  If  the 
reference  configuration  of  <R  is  specified  to  be  the  initial  configuration, 
say  at  time  t  =  0,  then  the  vector-valued  functions  on  the  right-hand  sides 

A  A 

of  (9.8)  can  be  identified  with  r(C,0)  and  d  (£,()),  respectively. 

'  » -  ■'*  -a 

Analogously  to  (9,6),  we  set 


If.  =  A  ,  If  -  (I)  , A  ) 
.■>  3  1  u  3 


(9.10) 


so  that  the  dual  of  (9.7)  is  given  by 


9 . 2  Basic  principles  of  a  Cosserat  curve  iR 

Consider  an  arbitrary  part  of  the  material  curve  S  in  the  present  con¬ 
figuration,  i.e.,  a  part  of  the  space  curve  l  bounded  by  f,  =  F  and 
f  =  f. 0  (F1  <  Fo)  »  ani1  let 

ds  =  (a33)V  ,  a33  =  a3.a,  (9.12) 

be  the  element  of  the  arc  length  of  SL .  It  is  convenient  at  this  point  to 

define  the  following  additional  quantities:  The  mass  density  p=p(f,t)  of 

* 

the  curve  £ ;  the  contact  force  n  =  n(F,t)  and  the  contact  director  forces 

ma  =  ma (f , t) ,  each  a  3-dimensional  vector  field  in  the  present  configuration; 

ct  ct 

the  assigned  force  f=f(F,t)  and  the  assigned  director  forces  £  =  £  (F,t), 
each  a  3-dimensional  vector  field  and  each  per  unit  mass  of  the  curve  £ ;  the 
intrinsic  (curve)  director  forces  k  =k  (F,t)  per  unit  length  of  £  which  make 
no  contribution  to  the  supply  of  moment  of  momentum;  the  inertia  coefficients 
yu  =  va(f)  and  yu^  =  ya^(f),  with  ya^  being  components  of  a  symmetric  tensor, 
which  are  independent  of  time;  the  specific  internal  energy  e=e(f,t);  the 
specific  heat  supply  r-r(f.,t)  per  unit  time;  and  the  heat  flux  h  =  h(F,t) 
along  J, ,  in  the  direction  of  increasing  f, ,  per  unit  time,  'Ifie  assigned  field 
f  represents  the  combined  effect  of  (i)  the  stress  vector  on  the  lateral 
surface  (2.20)  of  the  rod- like  body  denoted  by  f  ,  and  ( i i )  an  integrated 
contribution  arising  from  the  3-dimcnsional  body  force  denoted  by  f^,  e.g., 
that  due  to  gravity.  A  parallel  statement  holds  for  the  assigned  fields  in . 
Similarly,  the  assigned  heat  supply  r  represents  the  combined  effect  of 
(i)  heat  supply  entering  the  lateral  surface  (2.26)  of  the  rod-like  body  from 


the  surrounding  env i ronment,  denoted  by  r  ,  and  (ii)  an  integrated  contribution 

*  '  X 

The  notations  for  the  contact  force  n,  the  contact  director  forces  m  and  the 
curve  director  forces  k<l  differ  from  those  in  Green  and  Laws  (1966),  Green, 
Nagluli  and  Wenner  ( 1974a, b),  Naghdi  (1979a,  b  I  and  most  of  the  previous  papers 
on  the  subject.  In  fact,  the  vector  fields  n,ma,ka  of  Part  B  of  this  paper 
correspond,  respectively,  to  n,pn,7[a  of  Green,  Naghdi  and  Wenner  ( 1971a, b). 
.ifi!  most  of  the  previous  papers  on  the  subject 


arising  from  the  3-dimensional  heat  supply  denoted  he  r  .  Tims,  we  may  write 

,.  -  „a  .a  „a 

f  =  f  +  f  ,  l  =  9.,  *  i  ,  r=r+r 
~b  -c  -  ~b  ~-c  b  c 

We  assume  that  the  kinetic  energy  of  a  Cosserat  curve  $  per  unit 
length  of  the  curve  SL  in  the  present  configuration  is  given  by 

i  ,  ,  a  a8  , 

k  =  2P  v  •  v  ♦  2y  v  *  w  +  v  w  •  w, 

1  _  &  1  _  -u  -  -8 

We  further  define  the  momentum  corresponding  to  the  velocity  v  and  the 
director  momentum  corresponding  to  the  director  velocities  w  bv 

3k  ,  ct  ,  3k  ,  a  uP 

3V  =  P(!  +  y  ~a  *  3w“  =  tj(y  ^'  +  -v  W6J 

~a 

per  unit  length  of  Also,  the  physical  dimensions  of  p,n,f  arc 


phys.  dim.  p  =  (ML 


-2  -  > 
phys.  dim.  n  =  [MLT  ]  ,  phys.  dim.  f  =  |l,T 


where  as  in  section  3  the  symbols  |L],  [ M]  and  ( T ]  stand  for  the  physical 

dimensions  of  length,  mass  and  time.  The  dimensions  of  the  vector  fields 

m  and  k  depend  upon  the  physical  dimensions  of  d^.  Here  we  choose  d( 

Ct  ct 

to  have  the  dimension  of  length.  Then,  m  ,1  will  have  the  same  physical 

Ct 

dimensions  as  n,f  in  (y.16)  while  k  will  have  the  physical  dimension  of 

|MI.'V"  |  . 

Using  the  above  definitions  of  the  various  field  quantities  and  th 
not  at  ion 


|f(f  ,t)  L  =  f<f.,,t)  -  f If  .  ,t  )  , 


with  reference  to  the  present  configuration  the  conservation  laws  for  a 
* 

Cosserat  curve  are 


I  he  conservation  laws  (O.IX)  correspond  to  fqs.  (b.lo)  in  the  paper  of 
(.reen,  Nagluli  and  Wenner  (IP?4b). 


1  i) 


i  y .  i  ■» ) 


i  y .  i  r>  i 


(y.  u>) 


(9.  IS 


dt 


/  “p(v  +  yaw  )ds  =  /  "p  f  ds  +  [n | 

^1  '  M  '  ~  f'l 


"Jt  /jr^<A  +  yU6^6)ds  =  /^(P*“  -  (a33r  V)ds  +  [n.a]f'2  , 


,v  l,‘~p[r  x  v  +  y  (r  *  w  +  d  x  v)  +  d  x  y  w„]ds 
dt  ;  f  1  ~  y  ~  ~a  J  ~a  1  ~R 1 


1 ,  ■y 

=  / .  p(r  x  f  ■*•  d  x  £a)ds  +  [r  x  n  +  d  x  nia j 


a,  ^2 


-a 


u  "  '  F 


'1 


The  first  of  (9.18)  is  a  statement  of  the  conservation  of  mass,  the  second 
is  tlie  conservation  of  linear  momentum,  the  third  that  of  the  director 
momentum  and  the  fourth  is  the  conservation  of  moment  of  momentum.  We 
also  record  the  law  of  conservation  of  energy  in  the  form 

d  K  9  £  ->  (  f  , 

dt  /  “p  (c  +  Kjds  =  /_  “p  (r  ♦  f «  v  +  £a  •  w  )ds  +  [n  •  v  +  ma  •  w  -  h]  ~ 

*'l  9,  ~  -  ~  ~a  ~  '  ~  ~a  ■) 

The  basic  structure  of  (9.18)  and  (9.19)  are  analogous  to  the  cor- 

*  . 

responding  conservation  laws  of  the  classical  3- dimensional 

theory.  The  structure  of  (9.18)  and  (9.18).  are  less  obvious,  but  a 

j  4 

motivation  for  their  forms  is  provided  by  a  derivation  of  the  basic  field 

equations  for  rod-like  bodies  obtained  from  the  3-dimensional  equa- 

★ 

tions  of  continuum  mechanics  in  which  the  position  vector  r  in  3-space  is 
approximated  by  an  expression  of  the  form  (2.30).  It  should  be  noted  that 
the  conservation  laws  (9.18)  and  (9.19)  are  consistent  with  the 
invariance  conditions  under  superposed  rigid  body  motions,  which  ordinarily 
have  wide  acceptance  in  continuum  mechanics.  Moreover,  the  conservation  laws 
(9.18)  ,  (9.18),  and  (9.18)^  are  equivalent  to,  and  can  be  derived  from  the 
conservation  of  energy  (9.19)  and  the  invariance  conditions  under  superposed 


(9.19) 


rigid  body  motions.  The  conservation  law  (9.18)  for  the  director  momentum 
must  be  postulated  separately. 

Returning  to  the  conservation  laws,  after  making  suitable  continuity 
assumptions,  by  usual  procedures  from  (9. 18)  and  (9.19)  follow  the 

local  field  equations 

A  =  A(fJ  =  0(333)^  or  Pa33  +  pa3  •  =  0  ,  (‘d.JOJ 


Dn 

if  +  Xl  =  A^*y°wa)  , 


3  m 

TT 


.  .a  , a  ,  ,  a"  a8‘ 

+  A£  =  k  + A (y  v  +  y  w. 


and 


where 


Dd 

a,  x  n  +  d  »ka  +  x  ma  =  0 
~a  ~  at. 


Ar  -  ||  -  Ae  +  P  =  0  , 


3v  9w 

.  a  a  ~a 

P  =  n  •  -rzr  +  k  •  w  +  m  •  -3— 

dt,  ~  ~a  ~  9f 


is  the  mechanical  power. 


(9.21  ) 


(9.22) 


(9.27,) 


( 9.  21  ) 


(9.  2f>) 


9.7,  llie rare!)  i  cal  theories  of  Cosscrat  curves 

Although  the  theory  outlined  in  subsection  9.2  is  sufficiently  general 
for  many  applications,  on  occasion  it  becomes  necessary  to  consider  a  more 
general  theory  of  Cosscrat  curves.  Therefore,  we  now  briefly  discuss  the 
kinematics  and  the  balance  laws  of  Cosserat  curves  5?  having  I.  (  >  2)  directors 
attached  to  every  point  of  a  material  line  J,  the  number  I.  being  given  by 

r>.  1 ) . 


(.0, 


h 


is,  instead  of  (9.2)j  7,  we  specify  a  motion  offR^.  hv 


r  =  r(f,tj  d  =  d  lf  ,t)  (N=  1.- . K)  , 

where  the  vector  functions  d  are  assumed  to  he  symmetric  in  the 

~0lia2‘  '  ‘aN 

indices  a,av  .  .a...  The  velocity  vector  is  still  given  hy  (9.4)  but  cor 
1  Z  N  1 

responding  to  (9.4)7  we  now  define  the  director  velocities 


(  9  .  .’<!  ) 


~aiar--%  ~VY"hn 

We  recall  that  for  K=1  (#  =  A) ,  the  kinetical  quantities  and  the  assigned 

r  l  ■  i  1  It 

fields  introduced  in  subsection  9.2  consist  of  n,k  ,m  arid  f,£  -  Keeping  this 
in  mind,  for  a  body  we  admit  the  more  general  kinetical  quantities  and 


19.27) 


assigned  fields 


ala2 '  ‘  -0tN  ala2-'-UN 
n  ,  k  »  tn 


aia2 ' ' 'aN 


(  9 . 2  S  ) 


f  ,  a  , 

and  corresponding  to  (9.14)  and  (9.15)j  2  write  the  more  general  expressions 

for  kinetic  energy  of  and  associated  momentum  and  director  momenta, 

K 


name  1 v 


K  =  p  ( v  •  v  +  2  E  y 

~  ~  N=  1 


K  or^.-.o^ 


'•ala2-"UN 


K  ar..aNgr..eM 


+  E  y 

N=1 , M= 1 


w  •  w  | 

VY  br--RM 


(9.2'M 


K  a,... a 

4k  .  r  1  N  ! 

to  =  Plv  +  E  y  . • -a  1 

9v  ~  N=1  N 


ot.  . .  .a.,  a, .  .  f  •  ■  •  !\, 

------  -  p I y  Nv  •  >:  y  1  N  1  X  „ 

•*  M 


'*!  •  •  -(‘N 

each  per  unit  length  of  the  curve  l.  The  inertia  coefficients  v 

,V-Vi-liH.  1M  f  ,  'V  •  •% 

v  in  are  functions  of  f,  only,  y  art*  symmetric 

a  . .  .ai  .  .  .p  H ,  •  .  . B  u  .  -  .-t 

...  _  .  , .  1  N  1  M  1  Ml  \  , 

with  respect  to  indices  v  =y  and  are 

IN' 

also  symmetric  with  respect  to  UjU,. . .o^  and  fl  g,,. -B  .  In  the  special  case 
of  k-  1  ( fR  -iR|  ,  1 .  =  J ,  we  may  use  tiie  notations 


d  =  d 

-Uj  ~a 


w  =  w 

•Uj  -a 


a ,  B  . 
a  1  1  aS 

v  ,  v  =  v 


( t) .  s() ) 


for  a  detailed  statement  of  conservation  laws  appropriate  for  Cosserat 


curves  iR  we  refer  the  reader  to  Naghdi  ( 1 ;  >  7*  U  > ,  Sec.  .’ )  ,  hut  indicate  below 
K 


the  structure  of  the  corresponding  local  field  equations.  Thus,  for  the 
purely  mechanical  theory  of  Cosserat  corves  $  the  local  field  equations 


.  .  .a. 


=  0  ( N  =  1,...,2K)  ,  X  =  0  ,  A  =  XU)  -  -  i  a,. 


+  Xf  =  0  , 


a  ..,x  n  ...a 
.1  =  k  (  N  =  1 _ k  I 


K  >.x  .  .  .  .  ( i  1  x «  •  «  •  ■  < . .  ^  i  ■  •  ■  k 

,  ,  1  N  IN  I  % 

a .  '  n  ♦  \  Id  U  + - -- - 'in  i  o 

N-l  -‘l - *N 


whe  re 


/W  "aN  .  Yld' '  '‘n  V.1 . \ 

1  V  Mb 


((».  Pi 


■4 


k  ii ,  .  .  .  a 


t"  =  f  -  v  i'  y 
Ni=l 


w  (N  I . K )  , 

-'VaN 


li •  •  ar--uN-  J  ai---ltNpr--V 

=  y.  -  y  v  -  l  y  w 


M=  1 


V-A, 


( N  =  1  . . . . k| 


Also,  for  Cosserat  curves  the  expression  for  median  i  c:i  1  power  cor- 

K 


ri-^poikl  i  nj*  to  (^.2M  is 

*  v  k  .  .  .1 
h  V 

\=  I 


’)w 


n  •  .>  +  y  k 

if 


N 


k  V..«N  aN 

W  +  L.  HI  - . -rp - 

■^Ct  ...  I!  <)F 

1  N  N=  1 


Ihe  general  development  for  Cosserat  curves  91  outlined  ahove  is 
contained  m  a  paper  hv  Najjhdi  (1979,  Sec.  21.  When  k  =  1,  the  results  in 
s uh sect  ion  9.5  reduce  to  those  of  subsection  9.2  for  a  Cosserat  curve  iP . 


1. 1 . 


In.  1  I  a  s  t  u'  rod 


Within  tin.'  scope  of  the  theory  of  a  Cosserat  curve  iR  out  I  i  nod 
section  se  discuss  briefly  tile  constitutive  equations  for  elasti 
rod  in  the  presence  ot  tinite  detormat ion.  As  in  section  1,  we 
suppose  the  existence  of  a  strain  energy  or  stored  ener^v  per  unit 
s  t  I  such  that  ef  is  equal  to  tile  mechanical  poser  I’  defined  l> 

C.  Jo)  ,  i  -e.  , 


In  the  development  >f  nonl.near  . ->n  titutiie  equation  ;  tor  ela 
rods,  i^e  assume  that  the  strain  emir  dvn.it.  .  at  each  material  p 
ot  iR  and  for  all  t  is  specified  l>>  a  response  function  shich  depend 
r,d  in  I'd.J)  and  their  partial  deiivative.  uith  respect  to  Hut 

t  In.  response  function  must  remain  unaltered  under  superposed  rigid 
trail  National  di.pl  acemen  t  ,  the  dependenee  on  r  must  he  exeluded. 
tlie  constitutive  assumption  for  the  strain  energy  densits  can  In-  wr 
,1  •- 
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.  i  .u 

sh.  i'..  ,i  superposed  prire  stands  for 

1  i  '  1  I  /  •  •  , 

and  so  .  1 1  .u  make  similar  constitutive  assumptions  for  u.l  '  ,  in  in  (n 
In  t  In  e  v  on  ■■  t  i  ‘  ii  t  i  ve  equations,  which  represent  the  meihanic.il  ta 
ot  the  m  d  i  i  itii  .  the  dependence  of  the  response  functions  on  t 'u  i 

"■  ■■  n>  •  t  ric.il  properties  of  a  reference  state  and  material  i  nhor-ovi  m 
i  i  nd  1 1  a  t  ol  through  t  lie  argument  X. 

eon  t  i  t  ut  i  vc  tlivots  ><f  el.e-tu  i”if.  h 
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\  .'.em-r.il  development  of 


on  the  assumption  of  the  type  (10.J)  is  contained  in  the  paper  of  (Ireen , 
Naghdi  and  Kenner  (107-lb).  In  the  rest  of  this  section,  we  limit  the 
discussion  to  an  elastic  rod  which  is  homogeneous  in  its  reference  con 
figuration  and  suppose  also  that  the  dependence  of  the  response  functions 
on  the  properties  of  the  reference  state  occurs  through  tin  values  of  the 
f i nemat i cu 1  variables  in  the  reference  state.  Then,  in  place  of  (l(I.J), 
we  have 


<;j  =  iMr'.d  ,d;  ;  K  7 , l) ,  I) 7 ) 

-  -U  „Ct  -  ~  i 


with  similar  assumptions  for  n,k  ,m  .  After  substituting  (1(1.  -1 )  into 
(I  '.ll,  by  usual  techniques  we  obtain  the  following  forms  for  tin- 
const  i  t  lit  i  ve  equat  i  oils  : 


=  A  iit-  ,  ka  -  A  M- 


along  with  the  restriction 


•  O'-  Tl  . 

-  1  ad 

which  is  obtained  from  the  conservation  of  moim  nt  •:  ei-.tum  and  which 

restricts  the  response  function  y . 

he  do  not  discuss  here  the  reduced  forms  of  the  .move  constitutive 
equations  resulting  from  invariance  requirements  under  superposed  rigid 
body  motions,  but  for  such  reduction  refer  the  reader  to  Preen,  Naghdi  and 
Kenner  l  l  p74l> )  .  .lust  as  with  the  equations  of  motion,  it  is  necessarx  in 
.mill  i  cat  ions  to  specific  problems  to  obtain  alternative  forms  of  the  above 
constitutive  equations  or  their  reduced  forms  in  terms  of  tensor  compi  uent s . 
Such  component  forms  may  he  expressed  with  respect  to  liases  a.,  or  J  ,  m 
cor  respiuid  i  iig  i>ases  in  th.  reference  i  onfigur.it  ion  .  Reduced  forms  of 
l  I  - 1 .  ii  l  are  discussed  in  (ireen  et  a  I  .  (  1  uVlb  ,  See.  -|_ 


1  lo.-l ) 


1 1  o.(i ) 


1  1  .  Ide  ntification  of  the  assigned  fields  and  the  inertia  co  e  f  fi  cj  en  t  s 

The  local  field  equations  in  the  mechanical  theory  of  a  Cosserat  curve 

!R  have  the  same  forms  as  those  that  can  he  derived  from  the  3-d i mens ional 

field  equations  (2.9)  0  _  by  suitable  integration  over  the  cross-sectional 

1  a 

area  of  the  rod-like  body  with  respect  to  0  and  0“  [recall  the  definition  of 
a  rod-like  body  at  the  end  of  section  2]  and  in  terms  of  certain  definitions 
for  integrated  mass  density  and  resultants  of  stress  (for  details,  see  Green , 
Naghdi  and  Wenner  1974a).  Similarly,  the  energy  equation  (9.24)  has  the 
same  form  as  that  which  can  be  derived  from  the  energy  equation  in  the  3- 
dimcnsionul  theory  by  suitable  integration  over  the  cross-section  area  of 


the  rod-like  body  with  respect  to  0  and  0“  and  in  terms  of  certain  definitions 
for  integrated  internal  energy  density  and  heat  flux  in  the  3-dimensional 
theory  (see  (ircen  and  Naghdi  197n).  To  elaborate  further,  we  confine  atten¬ 
tion  to  t he  purely  mechanical  theory  and  recall  the  definitions 


A  =  pa  ‘  =  /  A  do'dQ"  ,  A  =  p  g^  , 
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and  the  expressions 
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,  I* , f  which  occur  in  ( 1 1  .  1  I  - ( 1  I . 4 )  are  defined  in  section  2  [lol- 

(2.9)  |,  the  line  integrals  are  taken  along  the  curve  f  =  const,  on 

orial  surface  (2. Jo).  A™ =  A  •  ga  and  A  -  \ftg  +  g^  is  a  vector  tangential 

*  'a  *  * 

surface  i  2 . 2o  I  so  that  A  •  o  =  A  v  +u,  =  (). 

rt  a 

wi  now  adopt  the  approximation  (2.  9ij,  then  there  is  a  I  1 


(  1  I  .  1  ) 
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correspondence  between  the  1 -dimensional  field  equations  that  follow  from 
the  conservation  laws  of  a  Cosserat  curve  and  those  that  can  be  derived 

from  the  3-dimensional  equations  provided  we  identify  r  and  the  director  d 

~<z 

in  (2.30)  with  (9.2)^  and  f  S) .  2  J  ^  >  respectively,  and  adopt  the  definitions 
(11.  M  (11.1),  as  well  as  the  definitions  of  the  resultants  mentioned  above . 

A  similar  1-1  correspondence  can  be  shown  to  hold  between  the  1  - d i mens i ona 1 
energy  equation  in  the  theory  of  a  Cosserat  curve  and  an  integrated  energy 
equation  derived  from  the  3-dimensional  energy  equation. 

The  various  quantities  in  (9.13)  are  free  to  be  specified  in  a  manner 
which  depends  on  the  particular  application  in  mind.  Also,  we  remark  that  in 
the  context  of  the  theory  of  a  Cosserat  curve,  the  inertia  coefficients 
y  , y  '  and  the  mass  density  p  require  constitutive  equations.  Indeed,  f  ,  8. 
and  r  ,  as  well  as  f^,£^  and  r^,  can  be  identified  with  corresponding  expres¬ 
sions  in  a  derivation  from  the  3-dimens.ional  equations  indicated  above  (for 
details,  see  Naghdi  19  9a).  Likewise,  p  and  the  coefficients  y°t,y°t^may  be 
identified  with  easily  accessible  results  from  the  3-dimensional  theory. 

In  what  follows,  we  assume  that  the  above  identifications  have  been 

ct  aB  ct 

made  and  that  the  quantities  p,y  ,y  ,  f  ,  are  known  or  specified.  The 
knowledge  of  f  ,  8.'^  depends  on  the  nature  of  the  boundary  conditions  on  the 
lateral  surface  of  the  particular  rod-like  body  under  consideration:  they 
may  be  specified  as  known  quantities  on  the  surface  (2.2b)  or  they  are 
unknown  and  must  be  determined  as  part  of  the  solution  of  the  problem. 


"1 


1 _ .  Additional  remarks  on  rods 

Topics  corresponding  to  those  in  Secs.  (>  and  7  have  so  far  received 
less  attention  in  the  case  of  rods  and  consequently  the  discussions  that  follow 
are  somewhat  brief.  We  first  consider  a  class  of  constraints,  apply  the  results 
to  an  incompressible  Cosserat  curve  and  then  go  on  to  briefly  comment  on  some 
recent  researches  which  bear  on  various  aspects  of  elastic  rods. 

Consider  a  class  of  constraints  which  are  linear  relations  between  the 
kinematic  variables 
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(N=  1,2,...,K) 


Similar  to  the  development  in  Sec.  (>  for  shells,  we  consider  (<}+!) 

J. 

constraint  equations  of  the  form 
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(12.1) 


k  Ma.av..a, 

+  Z  C  '  •  w/  =  0  (M=  o,  1 ,2 . Q)  , 

N=1  12  N 

where  A  ,11  ,C  “  are  vector  functions  of  d.  ,d.  only  and  do 

not  depend  explicitly  on  the  variables  (12.1).  We  assume  that  each  of  the 
nla2' ‘ -aN  a]a2< • *aN 

Inactions  n,k  ,m  arc  determined  to  within  an  additive 


(12.2) 


constraint  response  so  that 


n  =  n  +  n  , 
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are  specified  ny  consitutive  equations  and 
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(12.1) 


which  rep  resent  the  response  due  to  constraints  are  arbitrary  functions  nt 

I  lie  development  between  (  1  2 . 2  )  -  (  1 .’ .  (-> )  is  similar  to  that  for  mechanical 
,onsti  .in’s  in  the  .^-dimensional  theory. 
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t  and  are  workless.  Thus,  recall  mg  the  expression  ( () .  34 ) ,  we  set 
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(  13.!.) 


lor  all  values  of  the  variables  (12.1)  subject  to  the  constraint  con • 
c  1  i t ions  (12.2).  It  then  follows  that 
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where  p  =  p  (f, ,  t )  are  arbitrary  functions  which  play  the  role  of 
Lagrange  nil  1  tipi  i  ers  . 

We  consider  now  an  incompressible  Cosserat  curve  9  with  two  directors 
within  the  scope  of  the  above  constrained  theory.  As  in  the  case  of  an 
incompressible  shell-like  body  discussed  in  Sec.  6,  the  conditions 
representing  approximately  the  (3-dimensional)  incompressibility  condition 
(2.1(1)  may  be  derived  with  the  use  of  approximation  (2.27)  for  N=  1  given 
by  (2.30).  Under  this  approximation  (2.30),  the  base  vectors  are  given  In 
g  =  d  ,  g,  =  a  +0ad/,  where  a7  is  the  tangent  vector  to  the  curve  0U  =  0 
and  a  superposed  prime  is  defined  by  (10.3).  Then,  from  the  incompressibility 
condition  ( 2 . 1 0 ) j  we  obtain  an  approximate  expression  as  a  linear  function 
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o  f  "  , 


in  the  form 
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or'  equivalently  as 
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whore  in  (12.7)  and  (12.8)  use  is  made  of  the  notation  (9.6)  ,  and  (9.7). 

Wo  now  generate  three  conditions  representing  incompressibility:  One  of 

1  0 

those  is  obtained  from  integration  of  (12.8)  with  respect  to  ft  .O'-  over  the 

cross-section  of  the  rod-like  body  and  the  other  two  are  obtained  by  first 

multiplying  (12.8)  by  0^  (A  =  1,2)  and  then  integrating  the  resulting  eipia- 

1  2 

t ion  with  respect  to  0  ,ft  over  the  cross-section  of  the  rod-like  body. 
Those  three  conditions  can  he  written  as 


/  o  .it  fi.,./  ,.1,  ,u  ,  .  /  ,a.  u  ..•>  /  ,, 

v  +  !  y  d  +  y  [  (d  •  d  )d  -  (d^^  *  d  )d  ) )  •  w ^  +  y  d  •  w^  -  0  , 
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lor  an  incompressible  Cosserat  curve  under  discussion,  from  (12.2)  the 
constraint  conditions  arc 


/  „Ma  ,.Ma  /  ...  . 
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and  the  constrained  response  obtained  from  (12.6)  has  the  form 
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whore  p  ,p  ,p,  are  the  l.agrangc  multipliers.  Guided  hv  the  three  conditions 

.  .  ,  .M  ,,Mu 

(12.9)  and  (12.10)  for  A=  1,2,  we  select  the  vector- val ued  functions  A  ,1' 

(M,<  in  112.12)  and  (12.13)  to  have  the  special  values 
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Then ,  it  follows  from  (12.13)  and  (12.4)  that  the  expressions  for  tin- 
constraint  response  are  given  by 


The  applicability  of  the  theory  of  Cosserat  curves  is  not  limited  to 
only  elastic  rods  but  in  fact  can  be  applied  also  to  problems  of  fluid 
jets.  These  developments,  which  pertain  to  both  inviscid  and  viscous  jets, 
have  been  discussed  in  the  papers  of  Green  and  Laws  (11)68),  Green  (I97f>, 

I D7(> )  and  Naghdi  (1979b). 

A  constrained  theory  of  a  Cosserat  curve  with  two  directors  is  dis¬ 
cussed  by  Green  and  Laws  (1973)  and  includes  as  a  special  case  results  cor¬ 
responding  to  those  of  the  Rernoul  1  i -f.ulcr  beam  theory.  The  theory  of  small 
ile format  ion  superposed  on  a  large  deformation  of  an  elastic  Cosserat  curve, 
together  with  a  discussion  of  stability  problems  of  rods,  is  given  by  Green, 
Knops  and  Laws  (1908)  and  some  simpler  problems  in  the  context  of  the  nonlinear 
theory  of  rods  are  discussed  by  Lricksen  (1970). 

Hie  development  of  the  theory  of  Cosserat  curves  in  Sec.  9  is  carried 
out  within  the  scope  of  the  purely  mechanical  theory. 


In  earlier  work  on 


the  thermo-mechanical  theory  of  rods  by  direct  approach  (Green  and  Naghd 
1970),  only  one  temperature  field  was  admitted  and  this  allowed  for  the 
characterization  of  the  temperature  changes  along  some  reference  curve 
such  as  the  central  line  of  a  rod  in  the  (3-dimensional)  rod- like  body, 
but  not  for  temperature  changes  in  the  cross-section  of  the  rod.  The 
latter  effect  has  been  incorporated  recently  by  Green  and  Naghdi  (1979b) 
into  the  thermo-mechanical  theory  of  Cosserat  curves,  together  with 
appropriate  thermodynamical  restrictions  arising  from  the  second  law  of 
thermodynamics  for  rods. 


In  parallel  to  the  development  of  section  8  for  shells,  for  some 

purposes  it  is  convenient  to  have  available  the  basic  equations  of  a  Cosserat 

curve  in  a  direct  (coordinate- free)  notation  and  this  is  the  main  purpose  of 

the  present  section.  Just  as  in  the  case  of  shells,  we  shall  see  that  the  has 

equations  for  rods  in  coordinate-free  notation  are  very  similar  to  those  of 

the  corresponding  equations  in  the  3-dimensional  theory  and  thus  may  be 

more  suitable  in  the  discussion  of  general  theorems  or  in  the  developments 

which  parallel  those  in  the  3-dimensional  theory. 

We  introduce  the  notations  grad  and  Grad  to  denote  the  right  spatial 

and  material  gradient  operators,  respectively,  with  respect  to  the  position 

on  the  curve  c  in  the  current  configuration  and  on  the  curve  S  in  the 

K 

reference  configuration.  The  corresponding  divergence  operators  will  be 

denoted  by  die  and  Div,  respectively.  In  particular,  for  a  vector- va 1 ued 

* 

function  V(£,t),  we  write 


grad  V  =  V'®d3  ,  div  V  =  V'.  d3  , 
Grad  V  =  V'®D3  ,  Div  V  =  V'  •  D3  , 


( 


where  a  prime  denotes  partial  differ*,  ’tiation  with  respect  to  £  and  the 
symbol  0  denotes  tensor  product.  Also,  the  spatial  curve  gradient  operator 
is  defined  by 

grad^V  =  v'a3  ,  ( 

for  all  scalar- valued  functions  V(£,t). 

As  in  section  8,  we  introduce  a  measure  of  deformation  by  the  tensor  l; , 
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It  is  clear  that  the  notations  grad,  Grad,  div  and  Div  in  this  section 
stand  for  operators  with  respect  to  position  on  the  curve  c  and  need  not 
he  confused  with  the  similar  notations  in  section  8  for  surface  operators. 


name 1  v 


F  =  d .  ®  D1  =  Crad  r  +  da®  Da  ,  (13.3) 

and  in  view  of  the  notations  ((>.(>)  and  (9.10)  wo  observe  that 

FD„  =  d,  =  a,  ,  ID  =  d  (13.1) 

...3  .3  3  ~ot  ~a 

From  t he  definition  of  the  determinant  of  a  second  order  tensor  used  in 
section  X  (following  (8.4)]  and  the  conditions  (9.2),  and  (9.8),,  wo  obtain 


dot  F  =  |d  d7d^]/  (l)j  l>2D3J  >  «  •  (13.3) 

The  tensor  F,  a  linear  operator  on  vectors  in  3-space,  is  nonsingular;  and 
there  exists,  therefore,  the  inverse  deformation  gradient  F  *  defined  by 

!•  =  1).  09  d '  .  (  1  3,(.) 

,  j 

The  inverse  operator  F  1  transforms  vectors  in  the  present  configuration 
into  vectors  in  the  reference  configuration,  i.e., 

F“  1  d  =  1>.  (13.7) 


and  it  follows  that 

I-'1)-’  =  F  F'1  =  1  =  d.«d'  =  I).  <59  I)'  ,  (13.8) 

„  -.  -  -l  - 1 

where  1  is  the  unit  tensor  in  3- space.  We  also  introduce  here  the  gradient 
nl  tin  d  i  tec  t  o  rs  by 

C  -  Crad  d  =  d  '*>  I)3  .  I  13.9) 

1 1  .  (»  -St 

Koval  ling  the  definitions  (9.4)  .,  for  the  velocity  and  the  director 


74  . 


velocities,  as  well  as  (9.5),  we  have 


!•'  =  d.fcl)*'  =  d„®l)'>  +  d  &  l)a  =  v  '*  I,3  ♦  w  I)1' 
.  1  .1  -u  .  -  -u 


(i  =  d  '  6$  1)  =  w'®  I)J 

..a  -a  -»  ~a  ~ 


l;  l;  *  =  d,&d3  +  d  ®d°  -  grad  v  +  w  K>  dU 

~  ^  -  3  -  -a  ~  b  ..a 


*  1  'i, 

G  I-  =  w/?5d'  =  grad  w 

-cu  ex  &  -a 


The  formulas  ( 1 3. 3) - ( l 3 . 1 0)  represent  the  main  kinematieal  results  in 

terms  of  the  gradient  tensors  f',G  and  their  rates.  We  now  turn  to  kinetical 
h  ~  ~cx 

quantities  and  note  that  the  contact  force  n  and  the  contact  director  force 
m i ,  as  linear  functions  of  d3,  can  he  expressed  in  the  form 


n  =  d  *  N  d3  ,  ma  =  dl  kfV  , 

3.1  ~  ~  -  .i3  ~ 


;ith  the  second  order  tensors  N.M01  defined  by 


di,  N  =  n®dT  =  n  d.  «d, 

33  ~  -  ~3  -i  ...i 


df,  M06  =  mU  54  d  ,  =  m'ad.  <54  d.  , 
3.i  ..  ^.i  .  .i 


where 


i  ia  a  ,i 

n  =  n  •  d  ,  m  =  m  •  d 


Also,  it  is  convenient  to  introduce  a  tensor  K  through 


ka  =  d  *3  K  da 


dr’..  K  =  ku»d  =  k  lad .  fcd 

-  -  ~.n  -i  .u 


Vvhr  iv 


I 


k 


int 


(  13.  i<>) 


Before  proceeding  further,  we  recall  the  divergence  of  a  second  order 
tensor  field  T  is  tie  fined  by 

c  •  d  i  v  T  =  d  i  v  (  t'c  ) 


for  all  constant  vector  c,  where  superscript  T  denotes  transpose.  Applying 
the  above  definition  to  the  tensor  N  in  fJ3. 12)  j  ,  and  recalling  (‘).1J)|,  we 
obta  i  n 


d l V  ( N  e )  = d i v  j ( n  •  C )  d ,/d i _ I 
c  ••  -  C  —  -  .*>  5.* 


dll 


1  n  •  e)cL/d„‘„  I  '  •  a"'  =  •  c 

..  ,i.i  ...  ns 


with  a  similar  result  for  the  tensor  M  .  Thus,  we  have 


;>n 


;imrt 


d  i  v  M  -  ;  -  ,  d  i  v  M  =  -~'— 

c  ns  c  -  ns 


(  I  e> .  1  "  I 


With  the  use  of  (13.1),  the  kinematical  results  ( 1 3 . 8 ) - ( 1 3 .  9 ) ,  and  (I3.17| 
from  the  conservation  laws  (9.18)  follow  the  local  field  equations 


p  +  ,i  div  v  =  0 

d  i  v  N  +  p  f  =  p ( v  +  yaw  ) 

c-  '  ~a 

. .  ..ft  „a  .a  .  a'  ' 

d  l  v  M  +  p?.  -k  =  p(y  v  +  y  ) 

|N  +  k  +  M“(C  !•''  '  )  '  |  =  |N  +  K  +  Ma(C  f"  1  )  1  |  ' 


(13.18) 


As  in  corresponding  results  in  section  8,  the  last  statement  in  (13.18)  is 
similar  to  the  symmetry  of  the  stress  tensor  in  the  3-d i mens i ona 1  theory.  In 
particular,  it  may  be  observed  that  a  ,  *  n  ,  d  *  k  and  d  '  *  in*  are,  respectively, 
the  axial  vectors  of  a  ! |N  ■  N  *  |  .  a  3  _  |  K- K  1  |  and  a„‘ „  |  M  *  ( (1  l;  *)*  -  (<!  I  *  )(  M  * )  *  |  . 

1 1  i.i  .  -  .i.i  .  >t  " 


furthermore ,  in  terms  of  the  hinetical  quantities  N,M  ,K  in  (13.12)  and 
I  13.  IN)  aiul  the  rate  quantities  (13.10)  ,,  the  mechanical  power  becomes 

l  i  y  *- 

a^l’  =  tri IF1 (N+K)  +  ] ( I* ) ' }  . 

With  reference  to  constitutive  equations  for  elastic  rods,  instead  of 
the  kinematic  variables  used  in  section  10,  we  now  employ  the  variables 
(13.3)  and  (13.0).  Thus,  corresponding  to  the  constitutive  assumption 
I  I  "  •  I  )  ,  we  now  wr  i  te 


*  =  *(lM:a  ;  ’ 


win.  re 


(i  -  tirad  !i  =  li'Ol'f  , 

It-.u  .n  -it 


113.211 


along  with  similar  assumptions  for  N.k.Nf1.  Then,  with  the  use  of  (10.1  |,  (13.10; 
and  (13.20),  by  usual  technique  we  obtain  the  following  alternative  forms 
of  the  constitutive  equations: 


m“  --  0  (iv'<59d.) 

'Ill  •  /> 


3 . 2  2  I 


the  first  of  which  can  be  resolved  into 


(  I'  &  d_  I 


(13.  2  3  I 


\c  k  now  I  i’d  gmoii  t  .  Ilu-  ’Vault-  repot  t  id  in  re  were  obtained  in  the  course  of 
research  supported  In  tin  If'  .  t'tti,,  e  I  Naval  Kesearch  under  font  fact 
\  ( U  i ;  *  1  - 1  N  (  nils,  I'  re  1 1  c  t  \  1  •  ’< .  I  i  m  u  i  t  h  the  I  In  i  ve  rs  i  t  \  of  f  a  1  i  t  o  rn  i  a  , 

If  •  rk  e  I  e  v  . 
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